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Abstract 

For the integrable higher-spin XXX and XXZ spin chains we present multiple-integral 
representations for the correlation function of an arbitrary product of Hermitian ele- 
mentary matrices in the massless ground state. We give a formula expressing it by a 
single term of multiple integrals. In particular, we explicitly derive the emptiness forma- 
tion probability (EFP). We assume 2s-strings for the ground-state solution of the Bethe 
ansatz equations for the spin-s XXZ chain, and solve the integral equations for the spin-s 
Gaudin matrix. In terms of the XXZ coupling A we define £ by A = cos£, and put it 
in a region < £ < tt/2s of the gapless regime: —1 <A<1(0<£< it), where A = 1 
(£ = 0) corresponds to the antiferromagnetic point. We calculate the zero-temperature 
correlation functions by the algebraic Bethe ansatz, introducing the Hermitian elementary 
matrices in the massless regime, and taking advantage of the fusion construction of the 
i?-matrix of the higher-spin representations of the affine quantum group. 
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1 Introduction 



The correlation functions of the spin-1/2 XXZ spin chain have been studied extensively through 
the algebraic Bethe-ansatz during the last decade [H El [31 (U [51 [6] . The multiple-integral repre- 
sentations of the correlation functions for the infinite lattice at zero temperature first derived 
through the affine quantum-group symmetry [8] and also by solving the q-KZ equations [9|fTU] 
have been rederived and then generalized into those for the finite-size lattice under non-zero 
magnetic field. They are also extended into those at finite temperatures [TT]. Furthermore, the 
asymptotic expansion of a correlation function has been systematically discussed [12J. Thus, 
the exact study of the correlation functions of the XXZ spin chain should be not only very 
fruitful but also quite fundamental in the mathematical physics of integrable models. 

Recently, the correlation functions and form factors of the integrable higher-spin XXX spin 
chains and the form factors of the integrable higher-spin XXZ spin chains have been derived 
by the algebraic Bethe-ansatz method [T3"l [TH [To] . In the spin-1/2 XXZ chain the Hamiltonian 
under the periodic boundary conditions is given by 

1 L 

3=1 

Here cr" (a = X, Y, Z) are the Pauli matrices defined on the jth site and A denotes the XXZ 
coupling. We define parameter q by 

A = (q + q- 1 )/2. (1.2) 

We define r\ and ( by q = exprj and rj = i(, respectively. We thus have A = cos(. In the 
massless regime: — 1 < A < 1, we have < ( < tc for the spin-1/2 XXZ spin chain (11.11) . At 
A = 1 (i.e. q = 1), the Hamiltonian (11.11) corresponds to the antiferromagnetic Heisenberg 
(XXX) chain. The solvable higher-spin generalizations of the XXX and XXZ spin chains have 
been studied by the fusion method in several references [161 El HH1 US EDI EH E21 [23]. The 
spin-s XXZ Hamiltonian is derived from the spin-s fusion transfer matrix (see also section 2.6). 
For instance, the Hamiltonian of the integrable spin-1 XXX spin chain is given by 

^xxx = 2 ' ~ ' ^J'+ 1 ) 2 ) • (- 1 - 3 ) 

3=1 

Here Sj denotes the spin-1 spin-angular momentum operator acting on the jth site among the 
N s lattice sites of the spin-s chain. For the general spin-s case, the integrable spin-s XXX and 
XXZ Hamiltonians denoted T^xxx an< ^ ^xxz> respectively, can also be derived systematically. 

The correlation functions of integrable higher-spin XXX and XXZ spin chains are associ- 
ated with various topics of mathematical physics. For the integrable spin-1 XXZ spin chain 
correlation functions have been derived by the method of g-vertex operators through some 
novel results of the representation theory of the quantum algebras [211 [23 EH E?l [28] . They 
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should be closely related to the higher-spin solutions of the quantum Knizhnik-Zamolodchikov 
equations [TO]. For the fusion eight- vertex models, correlation functions have been discussed 
by an algebraic method [29]. Moreover, the partition function of the six- vertex model under 
domain wall boundary conditions have been extended into the higher-spin case [30J. 

In a massless region < £ < vr/2s, the low- lying excitation spectrum at zero temperature 
of the integrable spin-s XXZ chain should correspond to the level-fc SU(2) WZWN model with 
k = 2s. By assuming the string hypothesis it is conjectured that the ground state of the 
integrable spin-s XXX Hamiltonian is given by N s /2 sets of 2s-strings [31]. It has also been 
extended into the XXZ case [32J. The ground-state solution of 2s-strings is derived for the 
spin-s XXX chain through the zero-temperature limit of the thermal Bethe ansatz [T5]. The 
low-lying excitation spectrum is discussed in terms of spinons for the spin-s XXX and XXZ spin 
chains [3TL 132] . Numerically It was shown that the finite-size corrections to the ground-state 
energy of the integrable spin-s XXX chain are consistent with the conformal field theory (CFT) 
with c = 3s/ (s + 1) [531 EH E£ 36J. Here c denotes the central charge of the CFT. It is also the 
case with the integrable spin-s XXZ chain in the region < £ < 7r/2s [37J ESI EH]- The results 
are consistent with the conjecture that the ground state of the integrable spin-s XXZ chain with 
< C < vr/2s is given by N s /2 sets of 2s-strings j22E2E3E3E3HfllHDI32!. Furthermore, 
it was shown analytically that the low-lying excitation spectrum of the integrable spin-s XXZ 
chain in the region < ( < n /2s is consistent with the CFT of c = 3s/(s + 1) [4TJ H2]- In fact, 
the low-lying excitation spectrum of spinons for the spin-s XXX chain is described in terms of 
the level- A; SU(2) WZWN model with k = 2s g3]. 

In the paper we calculate zero-temperature correlation functions for the integrable higher- 
spin XXZ spin chains by the algebraic Bethe-ansatz method. For a given product of elementary 
matrices we present the multiple-integral representations of the correlation function in the re- 
gion < ( < it /2s of the massless regime near the antiferromagnetic point (( = 0). For an 
illustration, we derive the multiple-integral representations of the emptiness formation proba- 
bility (EFP) of the spin-s XXZ spin chain, explicitly. Here the spin s is given by an arbitrary 
positive integer or half-integer. Assuming the conjecture that the ground-state solution of the 
Bethe ansatz equations is given by 2s-strings for the regime of (, we derive the spin-s EFP for 
a finite chain and then take the thermodynamic limit. We solve the integral equations asso- 
ciated with the spin-s Gaudin matrix for < ( < tt/2s, and express the diagonal elements in 
terms of the density of strings. Here we remark that the integral equations associated with the 
spin-s Gaudin matrix have not been explicitly solved, yet, even for the case of the integrable 
higher-spin XXX spin chains [13] . We also calculate the spin-s EFP for the homogeneous chain 
where all inhomogeneous parameters £ p are given by zero. Here we shall introduce inhomoge- 
neous parameters £ p for p = 1, 2, . . . , N s , in §2.4. Furthermore, we take advantage of the fusion 
construction of the spin-s i?-matrix in the algebraic Bethe-ansatz derivation of the correlation 
functions [15]. 

Given the spin-s XXZ spin chain on the N s lattice sites, we define L by L = 2sN s and 
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consider the spin- 1/2 XXZ spin chain on the L sites with inhomogeneous parameters Wj for 
j = 1, 2, . . . , L. In the fusion method we express any given spin-s local operator as a sum of 
products of operator- valued elements of the spin-1/2 monodromy matrix in the limit of sending 
inhomogeneous parameters Wj to sets of complete 2s-strings as shown in Ref. [15]. Here, we 
apply the spin-1/2 formula of the quantum inverse scattering problem [JJ, which is valid at 
least for generic inhomogeneous parameters. Therefore, sending inhomogeneous parameters Wj 
into complete 2s-strings, we can evaluate the vacuum expectation values or the form factors of 
spin-s local operators which are expressed in terms of the spin-1/2 monodromy matrix elements 
with generic inhomogeneous parameters Wj. Here, the rapidities of the ground state satisfy the 
Bethe ansatz equations with inhomogeneous parameters Wj. We assume in the paper that the 
Bethe roots are continuous with respect to inhomogeneous parameters Wj, in particular, in the 
limit of sending Wj to complete 2s-strings. 

We can construct higher-spin transfer matrices by the fusion method [22, 23]. Here we recall 
that the spin-1/2 XXZ Hamiltonian f 1 1.11) is derived from the logarithmic derivative of the row- 
to- row transfer matrix of the six- vertex model. We call it the spin-1/2 transfer matrix and 
denote it by t (1,1) (A). Let us express by V^> an (£ + l)-dimensional vector space. We denote 
by T^' 2s \\) the spin-£/2 monodromy matrix acting on the tensor product of the auxiliary 
space V^> and the iV s th tensor product of the quantum spaces, [y( 2s )^ N ^ \y e call it of 
type (£, (2s)® Ns ), which we express (£, 2s) in the superscript. Taking the trace of the spin- 
£/2 monodromy matrix T^ e ' 2s \X) over the auxiliary space we define the spin-£/2 transfer 
matrix, t^ ,2s \X). For I = 2s, we have the spin-s transfer matrix t^ 2s ' 2s \\), and we derive the 
spin-s Hamiltonian from its logarithmic derivative. 

We construct the ground state ^f^) of the spin-s XXZ Hamiltonian by the B operators 
of the 2-by-2 monodromy matrix r^' 2s ^(A). As shown by Babujan, the spin-s transfer matrix 
f(2s,2s) ^ commutes with the spin-1 / 2 transfer matrix £ ( 1>2s ) (A) due to the Yang-Baxter relations, 
and hence they have eigenvectors in common [18]. The ground state \ipg 2s ^} of the spin-s XXZ 
spin chain is originally an eigenvector of the spin-s transfer matrix t^ 2s ' 2s ^(A), and consequently 
it is also an eigenvector of the spin-1/2 transfer matrix ^'^(A). Therefore, the ground state 
spin-s XXZ spin chain can be constructed by applying the B operators of the 
2-by-2 monodromy matrix T^ 1,2s ^(A) to the vacuum. 

We can show that the fusion i?-matrix corresponds to the i?-matrix of the affine quantum 
group U q (sl2)- We recall that by the fusion method, we can construct the R- matrix acting on 
the tensor product V® <g> QSl EZl QHl QSl 1201 12H 1221 12^- We denote it by R^ 2s \ In the 
affine quantum group, the i?-matrix is defined as the intertwiner of the tensor product of two 
representations V and W [HJ Hi3 H6] . D ue to the conditions of the intertwiner the i?-matrix of 
the affine quantum group is determined uniquely up to a scalar factor [37], which we denote by 
Rv,w- Therefore, showing that the fusion R- matrix satisfies all the conditions of the intertwiner, 
we prove that the fusion i?-matrix coincides with the i?-matrix of the quantum group, Rv,w- 
Consequently, for £ = 2s the fusion i?-matrix, i?( 2s ' 2s )(A), becomes the permutation operator 
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when spectral parameter A is given by zero. This property of the i?-matrix plays a central role 
in the derivation of the integrable spin-s Hamiltonian. It is also fundamental in the inverse 
scattering problem in the spin-s case [48J. 

There are several relevant and interesting studies of the integrable spin-s XXZ spin chains. 
The expression of eigenvalues of the spin-s XXX transfer matrix t( 2s ' 2s )(A) was derived by Babu- 
jan [TTl [TBI EE] through the algebraic Bethe-ansatz method. It was also derived by solving the 
series of functional relations among the spin-s transfer matrices [22]. The functional relations 
are systematically generalized to the T systems [19]. Recently, the algebraic Bethe ansatz 
for the spin-s XXZ transfer matrix has been thoroughly reviewed and reconstructed from the 
viewpoint of the algebraic Bethe ansatz of the ?7(l)-invariant integrable model [oTJl I5T"] . Quite 
interestingly, it has also been applied to construct the invariant subspaces associated with the 
Ising-like spectra of the superintegrable chiral Potts model [52] . 

The content of the paper consists of the following. In section 2, we introduce the R- 
matrix for the spin-1/2 XXZ spin chain. We then introduce conjugate basis vectors in order to 
formulate Hermitian elementary matrices E m,n in the massless regime where \q\ = 1. We define 
the massless higher-spin monodromy matrices T^ ,2s )(A) in terms of the conjugate vectors, after 
reviewing the fusion construction of the massive higher-spin monodromy matrices T^' 2s )(A) 
and higher-spin XXZ transfer matrices, t^' 2s '(\), for £ = 1,2, . . ., as follows. We express the 
matrix elements of T^^-NsW ^ n terms of those of the spin-1/2 monodromy matrix 7o^, L (A). 
Here ^^.^(A) is defined on the tensor product of the two-dimensional auxiliary space 
and the Lth tensor product of the 2-dimensional quantum space, 

(y(i)^L Here we recall 

L = 2s x N s . In the fusion construction [15], monodromy matrix T^^-n M acting on the N s 
lattice sites is derived from monodromy matrix 7q ^' 2 „. £ (A) acting on the 2siV s lattice sites by 
setting inhomogeneous parameters Wj to N s sets of complete 2s-strings and by multiplying it by 
the N s th tensor product of projection operators which project (y^)® 28 to V^. In section 3, 
we explain the method for calculating the expectation values of given products of spin-s local 
operators. We express the local operators in terms of global operators with inhomogeneous 
parameters w < ^ ,s,e \ which are defined to be close to complete 2s-strings with small deviations of 
0(e), and evaluate the scalar products and the expectation values for the Bethe state with the 
same inhomogeneous parameters w^ 2s ' e \ Then, we obtain the expectation values, sending e to 0. 
Here we note that the projection operators introduced in the fusion construction commute with 
the matrix elements of monodromy matrix T^ 1 ' ^ of inhomogeneous parameters wf s ' ^ with 0(e) 
corrections. In section 4 we calculate the emptiness formation probability (EFP) for the spin-s 
XXZ spin chain for a large but finite chain, and then evaluate the matrix S which is introduced 
for expressing the EFP of an infinite spin-s XXZ chain in the massless regime with ( < it /2s. 
Here we solve explicitly the integral equations for the spin-s Gaudin matrix, expressing the 
2s-strings of the ground-state solution systematically in terms of the string centers. In section 
5, we present explicitly the multiple-integral representation of the spin-s EFP. We also derive 
it for the inhomogeneous chain where all the inhomogeneous parameters £ p are given by 0. 
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For an illustration, we calculate the multiple-integral representation of spin-1 EFP for m = 1, 
(E 2,2 ), explicitly. In the XXX limit, the value of (E 2,2 ) approaches 1/3, which is consistent 
with the XXX result of Ref. [13J . In section 6, we present the multiple-integral representations 
of the integrable spin-s XXZ correlation functions. We express the correlation function of an 
arbitrary product of elementary matrices by a single term of multiple integrals. For instance, 
we calculate the multiple-integral representation of the spin-1 ground-state expectation value, 
(E ' x ), explicitly, and show that it is consistent with the value of spin-1 EFP in section 5, i.e. 
we show (E 1,1 ) + 2(E 2,2 ) = 1. Finally in section 7, we give concluding remarks. 



2 Fusion transfer matrices 



2.1 i?-matrix and the monodromy matrix of type (1, 1® ) 

Let us introduce the i?-matrix of the XXZ spin chain [TJ, El SI E]- We denote by e a,b a unit 
matrix that has only one nonzero element equal to 1 at entry (a, b) where a, b = 0, 1. Let V\ 
and V2 be two-dimensional vector spaces. The i?-matrix acting on V\ ® V% is given by 



B+(A!-A 2 ) 



E 

a,b,c,d=0,l 



[U. 



cd 



b e a,c : 



J>,d 



( 1 



\0 





b{u) 
c+(u) 






c-(u) 
b(u) 




0\ 




1 / 



(2.1) 



where u = Ai — A2, b(u) = sinhu/ sinh(w + if) and c (u) = exp(±w) sinhr// sinh(u + rj). In the 
massless regime, we set 77 = i( by a real number (, and we have A = cos£. In the paper we 
mainly consider the region < ( < it /2s. In the massive regime, we assign r] a real nonzero 
number and we have A = cosh 77 > 1. Here we remark that the -R + (Ai — A2) is compatible with 
the homogeneous grading of U q (sl 2 ), which is explained in Appendix A [15] . 

We denote by R^ p \u) or simply by R(u) the symmetric i?-matrix where 0^(11) of (12.1ft are 
replaced by c(u) = sinhr// sinh(u + 77) [T5J. The symmetric i?-matrix is compatible with the 
affine quantum group U q (sl2) of the principal grading [15] . 

Let s be an integer or a half-integer. We shall mainly consider the tensor product V^ 2 ^ <8> 
• • '®Vn^ °f (2s + l)-dimensional vector spaces vj 2 ^ with L = 2sN s . In general, we consider the 

tensor product Vq 2s °^ (g> p^ 2si ^ ® • • • <g> Vr 2Sr ^ with 2si H h 2s r = L, where have spectral 

parameters Aj for j = 1, 2, . . . , r. For a given set of matrix elements A^'p for a, 
and a, j3 = 0,1, ... , 2s k, we define operator Aj k by 



0,1, 



2si 



Ajk 



a,b=l a,P 



(2s ) 



(2si) ® • • • ® /f_r° 



r(2a*_i) 
'fc-1 



a,/9 



r(2s fc+ i) 
J fc+1 



j(2Sr) _ 



(2.2) 



We now consider the (L + l)th tensor product of spin-1 /2 representations, which consists 
of the tensor product of auxiliary space Vq and the Lth tensor product of quantum spaces 
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y. (1) for j = 1,2,..., L, i.e. V (1) ® fy/^ ® • • • <g> vf^J . We call it the tensor product of type 
(1, and denote it by the following symbol: 



1,1^) = (1,1,1,...,!). (2.3) 




Applying definition (12.21) for matrix elements R(u)^ of a given i?-matrix, we define R- 
matrices Rjk(\j, \k) = Rjk(Xj — Xf.) for integers j and k with < j < k < L. For integers j, k 
and i with 0<j<k<£<L, the i?-matrices satisfy the Yang-Baxter equations 

Rjk(^j — ^k)Rji(^j ~ ^e)Rki(^k — Xg) = Rke(^k — ^i)Rje(^j — Xi)Rjk(Xj — Xf.) ■ (2.4) 
We define the monodromy matrix of type (1, 1® L ) associated with homogeneous grading by 
7 o,i2"'2(^o; w i, i» 2 ,..., w L ) = Rq L (\ -w l )--- Rq 2 (X - w 2 )R^(X - wi) . (2.5) 

Here we have set Xj = Wj for j = 1,2, ...,L, where Wj are arbitrary parameters. We call 
them inhomogeneous parameters. We have expressed the symbol of type (1, 1® L ) as (1, 1) in 
superscript. The symbol (1, 1+) denotes that it is consistent with homogeneous grading. We 
express operator- valued matrix elements of the monodromy matrix as follows. 

V C\ 2 ..! L {X] {wj} L ) D\ 2 ..' L (X;{wj} L ) J 

Here {wj}l denotes the set of L parameters, w±, w 2 , ■ ■ ■ , Wl- We also denote the matrix elements 
of the monodromy matrix by \T^-^t\{X\ {wj}i)}a,b for a, b = 0, 1. 

We derive the monodromy matrix consistent with principal grading, Tq^; l (X; {w^l), from 
that of homogeneous grading via similarity transformation Xoi—l as follows [15j. 

T oa2*l( X > { w j}l) = Xoi2- L T^;l p) L (X; {w j } L ) Xm 1 2 ... L 

X i2-lA[ 1 2 p) . l (X; {w 3 } l )xiLl e^^S^A; {w 3 } l )xi1.l 



e 



x »Xi2-lC£ p) l (X- { Wj } L )xll. L Xu-lD^ l (X; {w s } l )x&~i 



(2.7) 



Here Xoi -l = $o$i • • • and $j are given by diagonal two-by-two matrices $j = diag(l, exp(wj)) 
acting on Vj for j = 0, 1, . . . , L, and we set wq = Xq. In Ref. [15] operator A^ 1+ ^(X) has been 
written as A + (X). Hereafter we shall often abbreviate the symbols p in superscripts which 
shows the principal grading, and denote (2s p) simply by (2s). 

Let us introduce useful notation for expressing products of _R-matrices as follows. 

p(«0 _ td(w) t,(w) td(w) 

-'H^S-n — -"-Ira ■■'- rl 13 -"-12 i 

rl 12-- n-l,n — -R-ln "' rL n-ln- 

Here R^ denote the .R-matrix R^ = R^(X a — A&) for a, b = 1, 2, ... , n, where w = + and 
w = p in superscripts show the homogeneous and the principal grading, respectively. Then, 
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the monodromy matrix of type (1, 1® L w) is expressed as follows. 

For instance we have ^{^(Aq; {wj} L ) = [Rq^1. l (Xo] {wj} L )}o,i ■ 

2.2 Projection operators and the massive fusion R- matrices 

Let V\ and V 2 be (2s + l)-dimensional vector spaces. We define permutation operator 11^2 by 

III, 2 v 1 <g) v 2 — v 2 ® v\ , vx G Vx , v 2 G V 2 . (2.10) 

In the case of spin-1/2 representations, we define operator Rf 2 (\i _ A2) by 

i?+(A! - A 2 ) = Ux^RUXx - A 2 ) . (2.11) 

We now introduce projection operators P[ 2 ...t f° r ^ > 2. We define by = R\2{ r i)- 
For £ > 2 we define projection operators inductively with respect to £ as follows [4"B"j |2"5]. 

Pfh = pVIi^UM ~ 1 ■ (2-12) 

The projection operator P[ 2 ... t gives a g-analogue of the full symmetrizer of the Young operators 
for the Hecke algebra [46]. We shall show the idempotency: (P^...g) — Pi2--e ^ n Appendix B. 

Hereafter we denote -P 12 ...^ also by P{ for short. 

Applying projection operator Pa% 2 ... ai to vectors in the tensor product ® - ■ -®Val\ 

we can construct the {£ + 1) -dimensional vector space Va% 2 ... ai associated with the spin-£/2 
representation of U q (sl 2 ). For instance, we have PaiLl H — }a = (?/[2] ? )||2, l) a , where we have 
introduced | H — ) a = |0) ai (g) |l)a 2 - The symbols such as g-integers are defined in Appendix C. 
Moreover, the basis vectors \\£, n) (n = 0, 1, . . .,£) and their dual vectors (£, n\\ are given for 
arbitrary nonzero integers £ in Appendix C. We denote Vafl 2 ... ai also by 14 or for short. 

Since Pf 2 -i * s consistent with the spin-£/2 representation of U q (sl(2)) (see (IC.6P ). we have 

e 

^ = X>,n>&*ii- ( 2 - 13 ) 

n=0 

Applying projection operator ^ 2 §_ 1)+1 ...2 s(i _i) + 2 S to tensor product • • -®v£l_ 1)+2a > 

(2s) (2s) 

we construct the spin-s representation V 2s ^_ 1 - )+1 ... 2s ^_ 1 - )+2s . We denote it also by Vj , briefly. 
In the tensor product of quantum spaces ® • • • Cg) ■ , we define P^ S ] L by 

p(2») _TTp (2s) (2 14) 

12 - L — H^sCi-lJ+l • 
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Here we recall L = 2s N s . We have put 2s in place of £. 

We now introduce the massive fusion i?-matrix Rq'^ s + ^ on the tensor product Vq (g> vj 2s ^ 
(j = 1,2, . . . ,N S ). It is valid in the massive regime with A > 1. We first set rapidities A a 



of auxiliary spaces V^- by A Qfe = A ttl — (k — l)rj for k = 1,2, ...,£ — 1, and then rapidities 
A 2s (j-i)+fc of quantum spaces ^ 2 ( s 1 ( ) j _ 1)+fc by A 2s (i-i) + fc = A 2s (j-i)+i - (k - l)rj for k = 1, 2, . . . , 2s 
L,2,... 



and j = 1, 2, . . . , N s . We define the massive fusion i?-matrix Rq' 2s+S) as follows. 



R, 



l(\ _ \ \ — p{t) p( 2s ) E>+ pW p( 2s ) 

V^di ^2s(j-l)+i; - r a r --a £ - r 2s(j-l)+l X ai---a f , 2s(j-l)+l— 2sj - r ai---G^- r 2s(j-l)+l 

pW p( 2s ) E>+ . . . E>+ E>+ pW p( 2s ) 

J oi -o^- 1 2s(j'-l)+l oi-a«,2sj -"'ai---a ( .,2s(jr'-l)+2 L a 1 - -a t ,2s(j-l)+l ai-a* 2s(j-l)+l ' 



(2.15) 



2.3 Conjugate vectors and the massless fusion /^-matrices 

In order to construct Hermitian elementary matrices in the massless regime where \q\ = 1, we 
now introduce vectors \\£, n) which are Hermitian conjugate to (£, n\\ when \q\ = 1 for positive 
integers i with n — 0, 1, . . . , £. Setting the norm of | \£, n) such that (£, n\ | | \£, n) = 1, we have 



IK,") 



-(h+-+i„)+ni-n(n-l)/2 



l<h<--<i n <t 



' £ ' 


q -n(i-n) 








n 








1 





(2.16) 



Here we have denoted the binomial coefficients for integers £ and n with < n < I as follows. 

- I! 
// / {£ — n)\n\ 



(2.17) 



The g-binomial coefficients are defined in Appendix C. Dual vectors (£, n\\, which are conjugate 
to \ \£, n), are defined in Appendix C, and we have 

-2 



(£,n\\\\£,n) 



' £ ' 


2 




n 


g 


yn) 



(2.18) 



They are determined by the action of with opposite coproduct: A op = r o A. For instance, 
we have \\£, n) = const. A op (X~) n \\£,0) . Here and A op are defined in Appendix A. 

For an illustration, in the spin-1 case, the basis vectors ||2,n) (n = 0, 1,2) are given by [15] 

112,0) = | + +), 

||2,1) = \ + -) +q -i\-+), 

||2,2) = |--). (2.19) 

|1)2, briefly. The conjugate vectors ||2,n) (n = 0, 1,2) are given by 

Eoj = 1 + +), 



Here | H — ) denotes |0)] 



[2], 



2q ' 



|2,2) 



-) 



(2.20) 
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In the massless regime, operator ||2, 1)(2, 1|| is Hermitian while ||2, 1)(2, 1|| is not. 
Let us now introduce another set of projection operators -Fy...^ as follows. 

i 

n=0 

Projector Pf}^ is idempotent: (P^} e ) 2 = P^l t In the massless regime where \q\ = 1, it is 
Hermitian: (p^.g) = P\%- From (I2.13P and (I2.2ip . we show the following properties: 

r V2-l r \-l ~~ r \2-li \&.&&) 

P^Al, = (2-23) 
In the tensor product of quantum spaces, 1/^ (g) • • • ® V^ 2,s \ we define P^ S . L by 



4 2 i=n<ti) + i- ( 2 - 24 ) 

Here we recall L = 2sA^ s such as for (12.141) 
We define the massless fusion .R-matrix 
of conjugate vectors to the product of i?-matrices, as follows 



We define the massless fusion .R-matrix R^'^ s+ \ applying projection operators P consisting 



p(t> 2s +)(\ \ _ p(Q p(2s) p+ p(i) p(2s) 

J^Oj V^m UJ 2s (j-l)+l) — r a 1 -a e r 2s{j~l)+l IX a 1 -a l ,2s{j-l)+l-2sj r a 1 -a e r 2s{j-l)+l 

_ 5W p( 2s ) p+ .../?+ R+ pW p( 2s ) 

- 1 oi— at 2s(j-l)+l ±L a 1 --a l ,2sj 1 L ai ---at, 2s(j-l)+2 V ■■■a £ , 2s(j-l)+l tti— o/" 4 2s(j-l)+l ' 

(2.25) 

We should remark that the massless fusion .R-matrix _R( £ > 2s ) and the massive fusion .R-matrix 
jf[i,2s) jjayg the same matrix elements. Some examples are shown in Appendix D. 



2.4 Higher-spin monodromy matrix of type (£, (2s) s ) 

We now set the inhomogeneous parameters Wj for j = 1,2, ...,L, as iVs sets of complete 



2s-strings [151 . We define iu^I v , fl for /3 = 1, . . . , 2s, as follows. 



(2s) 
W 2s(6- 



i) +/ s = 6-(/3-l)»?, for 6 = 1,2,. (2.26) 



We shall define the monodromy matrix of type (1, (2s)® ") associated with homogeneous grad- 
ing. We first define the massless monodromy matrix by 



^(l,2s+) ( , r fiX s _ p(2s) .r,.,(2«)-| Np(2s) 

^ (2s+) (A;feU) f (2s+) (A;te}ivJ 
^ (2s+) (A;{a}/vJ -D (2s+) (A; {£&}jvj 



(2.27) 
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Here, the (0,0) element is given by A^iX; = P^ s) . L A^ + \X; {wf s) } L )p£ s } L . We then 

define the massive monodromy matrix by 

T (l,2s+)/ x \ p (2s) (1,1+),. , r (2s)-. n p (2s) 

^ (2s+) (A;{&K) B^+\\;{Z b } Ns ] 
C< 2s+ >(A; ^ (2s+) (A;{6W 



(2.28) 



Let us introduce a set of 2s-strings with small deviations from the set of complete 2s-strings. 

(2s; e) 
; 2s(6-l)+/3 



^(b-i) +0 = ^-(P-l)v + ^ ) , for 6=l,2,---,iV s , and (5 = 1, 2, . . . , 2s. (2.29) 



Here e is very small and rf^ are generic parameters. We express the elements of the monodromy 
matrix T^ 1 ' 1 ) with inhomogeneous parameters given by wf s,e ^ for j = 1, 2, . . . , L as follows. 

T , 12 ... L (\i^ M-[c£%<\\) D^(X))- (2 - 30) 
Here we recall that Af^l^iX) denotes A^} L (X\ {w^ s, ^}l)- We also remark the following. 

^2 S .1(A; feU) = KinP£°} L A^(X; {wf s ^} l )p£;} l . (2.31) 

Let us express the tensor product (g) (v} 2s ^ (g) • • • ® j , by the following symbol 



(€, (2a)®"-) = & 2a, 2a,..., 2a). (2.32) 

Here we recall that V^p abbreviates Vafl 2 ... av In the case of auxiliary space V^p we define the 
massless monodromy matrix of type (£, (2s) 0Afs ) by 

f(.t,2s+) _p W f(l,2s+) (x ^(1,2.+) ( , _„N...f(l,2s+) /■> _f/_iy 

J 0,12---Af s — ^ai^-a* 1 ai ,12---N s \ A ai J 1 a 2 ,12---Ns\ Aa i 'D 1 a e , 12-N S \ Aa l \ l 1 I'D r w aia 2 -a* j 

(2.33) 

and the massive monodromy matrix of type (£, (2s)® Ns ) by 

T (A2s+) _ pW T (l,2s+) fA x T (l,2s+) fA _ v (l,2s+) pW f2 34 x 

"'0, 12---AT S ~~ - r aia 2 ---a i . 1 a 1 ,12---N s \ A ai)- L a2 , 12---AT S '/J - 1 a e , 12-N S \ A ai \ C III a\a 2 ---ag " 

For instance, the (0, 1) element of the massive monodromy matrix T( 2,2s+ )(A) is given by 
(2,011^+^(^112, 1) = A%+\X)B% +)( A -r)) +q- 1 B£ + \\)Ag+\\ -v) . (2.35) 

2.5 Series of commuting higher-spin transfer matrices 

Suppose that \l,m) for m — 0, 1, . . . , I, are the orthonormal basis vectors of and their 
dual vectors are given by (£,m\ for m — 0, 1, . . . , £. We define the trace of operator A over the 
space V® by 

l 

\x v wA = ^2{£,m\A\e,m) . (2.36) 

m=0 
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The trace of A over V™ is equivalent to the trace of A over the £th tensor product of 
(yW)®^ multiplied by a projector P^' (or P^) as follows. 

ti vW A = tr (V (i ))9 i (P W A) 

= e (^):::- ( 2 ^) 

ai,...,a^=0,l 

It follows from (12.131) that the trace with respect to is given by (12.361) . 
We define the massive transfer matrix of type (£, (2s)® Ns ) by 

= E«^ n ll^;S.S.(A)2£S.S.(A - i?) • • -T^l(\ -{£- 1) V ) \\£,n) a , (2.38) 

n=0 

and the massless transfer matrix of type (£, (2s)® Ns ) by 
*S-Sv?(*) = tr v(') (^ofi?.Jr.(^)) 
= £a& n ' I f i^ 2 2i(A)T a ( 2 1 ;? 2 S l s (A - r/) • ■ ■ (A - (* - 1)^) jjl^) a . (2.39) 

n=0 

It follows from the Yang-Baxter equations that the higher-spin transfer matrices commute 
in the tensor product space Vy ® • • • ® , which is derived by applying projection operator 
Pi2 S ..L ^° ® ' ' ' ® ■ For instance, for the massless transfer matrices, making use of (12.221) 
and (I2.23P we show 

P£?M* S N + JW, = , for £, m G Z> . (2.40) 

Therefore, for the massless transfer matrices, the eigenvectors of ^-iv" W constructed by 
applying B^ 2s+ \\) to the vacuum |0) also diagonalize the higher-spin transfer matrices, in 
particular, t^'^N G\)- Thus, we construct the ground state of the higher-spin Hamiltonian in 
terms of operators B^ 2s+ \\), which are the (0, 1) element of the monodromy matrix T^' 2s+ \ 



2.6 The integrable higher-spin Hamiltonians 

We now discuss the integrable massless spin-s XXZ Hamiltonian. For (2s + l)-dimensional 
vector spaces Vj_ and V 2 \ we can show that the massive spin-s fusion .R-matrix R^' 2s + ^ (u) 
at u = becomes the permutation operator n 12 for ® V^ 2 ^ . Furthermore, operator 

Ri2 2s (u) = U.i j 2Ri2 2S+ \ u ) has the following spectral decomposition: 

R?T S+ \u) = f> s - 2 » Ut%).., , (2-41) 

3=0 
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where operator {P2l-2j)i,2 projects V} ®V 2 to spin-(2s— j) representation for j = 0,1, ... ,2s. 
Functions P4, s -2j(u) are given by [15] 



. . -^4- sinh(kri — u) , 

P4s- 2J « = n • h(h ' • 2 - 42 

_- LJ -. smn(A;r7 + u) 



The massless spin-s .R-matrix is thus given by 

CfrV) = £ Pis-M P£g_ 1)+1 Pffl a • (nM) M+1 • (2-43) 

It is easy to show that the massless spin-s i?-matrix i?[ 2 2 s,2s+ '(-u) becomes the permutation 
operator at u = 0: i?[ 2 2 ' 2s ^ (0) = n 1;2 - Therefore, putting inhomogeneous parameters £ p = 
for p = 1,2, ... , N s , we show that that the transfer matrix t[ 2s .'. 2 ^ s + ' ) ( A) becomes the shift operator 
at A = 0. We thus derive the massless spin-s XXZ Hamiltonian by the logarithmic derivative 
of the massless spin-s transfer matrix, similarly as for the massive case. 

N 3 



nj{2s) _ d , ~( 2s ,2 S +) m 



A=0,fc=0 i=x C/ " 



d ~(2 S ,2 S ) 



M=0 



N s 2s , 
i=l j=0 

3 Higher-spin expectation values 
3.1 Algebraic Bethe ansatz 

In terms of the vacuum vector |0) where all spins are up, we define functions a(A) and d(X) by 

A^\X;{ Wj } L )\0) = a(X;{ Wj } L )\0), 

D(^)(A;K} L )|0) = d(X;{ Wj } L )\0). (3.1) 
We have a(A; {wjji) = 1 and 

d(A;WL) = n^.«'j)- (3-2) 

3=1 

Here b(X, fi) = b(X — fi). For the homogeneous grading (w = +) and the principal grading 
(w = p), it is easy to show the following relations: 

A( 2s -)(A)|0) = ^ 2sw \X)\0} = a ( - 2s \X;{^})\0}, 

D {2sw) (X)\0) = D {2sw) (X)\0) = S 2s) (X;{£ k })\0) , (3.3) 
where a^ 2s ^(A; {6c}) and d^- 2s \X; {£&}) are given by 

« (2s) (A;fc}) = o(A; {wf s) }) = 1 , 

^ (2s) (A;fc}) = rf(A;{^f ) }) = n^(A,e P ). (3.4) 

p=i 
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Here we have defined b t (X, //) by b t (X, fi) = sinh(A — //)/sinh(A — p, + trj) . Here we recall 
b(u) = bi(u) = sinhw/ sinh(-u + 77). 

In the massless regime, we define the Bethe vectors |{A Q } M ) for w = + and p, and their 

dual vectors ({A a } M | for w = + and p, as follows. 

\{K)m W) ) = f[B {2sw \\ a )\0) , (3.5) 



-(2sw) 



a=l 

M 



I — LL ' / "1 — r ' — ' / \ 

({K} M I = (0\i[C^\\ a ). (3.6) 

Q=l 

Here we recall B^ 2s+ \\ a ) = P[ 2s) L B^ + \\ a ,{w k } L )P (2s) L . The Bethe vector (E3D gives an 
eigenvector of the massless transfer matrix 

^ 2sw \^ = A<*"»>Qi; + 5( 2 ^)(/i; {£,U) (3.7) 



for ?u = + and w = p with the following eigenvalue 

m . , x TV. M 



if rapidities {Xj}m satisfy the Bethe ansatz equations 

b(Xki A/ 



sinh(Aj — fi + rj) t yj u ( t ^ yj sinh(/i — Aj + rj) 



(3.* 



n^(^^) = nSrxy = AO- (3-9) 

p=l fc^j ^ -" ' 

Let us denote by |{A a (e)}^ s *"'^) the Bethe vector of M Bethe roots {Aj(e)}M for w = +,p: 

\{U^} { M W]e) ) = S (2awie) (Ai(e)) • • • 5 (2 ^ ;e) (AAf(e))|0) , (3.10) 

where rapidities {A^e)}^,/ satisfy the Bethe ansatz equations with inhomogeneous parameters 
wf s as follows. 

It gives an eigenvector of the transfer matrix 

t^ lw \p; {™f ;e) }i,) = A^'fa {«,f ;e) } L ) + D^Gu; ;e) } L ) (3.12) 
with the following eigenvalue: 

J a = a smh(A i (e) - p) ax sm h(/x - A^e)) 

(3.13) 

Let us assume that in the limit of e going to 0, the set of Bethe roots {\j{e)}M is given by 
{\j}u- Then, we have 

P§Sl \Mm +) ) = }™ PISl l{A,(e)}S s+;e) ) • (3.14) 
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3.2 Hermitian elementary matrices ^ TO ' n ( 2s+ ) [ n ^he massless regime 

We define massless elementary matrices E m ' n( - 2s+ ^ for m, n = 0, 1, . . . , 2s, in the spin-s repre- 
sentation of £7g(sZ 2 ) as follows. 

E m ' n( - 2s+ ) = \\£m)(i,n\\. (3.15) 
In the tensor product space, (y( 2s ^)® Ns , we define _g™-' n ( 2s +) for i = 1, 2, . . . , iV 8 by 

£m,n(2s+) _ ^-(2s)^(i-l) ^ _gm,n(2s+) ^ ^j(2s) j®(2V s -i) _ (3.16) 

Elementary matrices _E n > ri ( 2s +) for n = 0, 1, ... , 2s, are Hermitian in the massless regime. In 
fact, when \q\ = 1, for m, n = 0, 1, . . . , 2s, we have 

-l 



^JTj m , n (2s +) 



i 



2s 
m 



2 r -1-2 

2s 
n 



2s 
m 



2s 
n 



EJ't, (2s +) 



(3.17) 



We can express any given spin-s local operator of the massless case in terms of the spin-1/2 
global operators by a method similar to the massive case [15]. For m = n, we have 



E 



n, n (2s+) 



2s 
n 

x 



2s 

n 



(«-l)2s n 

p?i n (^ (i+) +^ (i+) )K)n^ (i+) K-^) 

a=l k=l 
2sN s 

w(i-D2 S+k ) n (^ (i+) +^ (i+) )K)p 1 (2 i 



(3.18) 



fc=n+l 



a=i2s+l 



Formulas expressing i?" 1 ' n ( 2s+ ) for m > n or m < n are given in Appendix E. 

When we evaluate expectation values, we want to remove the projection operators intro- 
duced in order to express the spin-s local operator in terms of spin-1/2 global operators such 
as in (I3.18p . Then, we shall make use of the following lemma. 

Lemma 3.1. Projection operators P^1 L and P^1 L commute with the matrix elements of the 
monodromy matrix 7oi2-^!(^i {wf S, ^}L) such as A^ 2s+ ' ,e \\) in the limit of e going to 0. 

P£Lrfd+kM*t' £) h) P ™L = ^i*2 1 ^(A;{^ ;e) }L) + 0(e), (3.19) 



P£l^+k*^' e) }L)Pg:} L = P[i! L T^...^A^f S ' e) } L ) + 0{e). (3.20) 
For instance we have P^} L B ( - 2s +' e \X)P^ s) . L = P^} L B^ 2s+ > e \\) + 0(e). 
Proof. Taking derivatives with respect to inhomogeneous parameters Wj, we can show 

zfe^iCA; H^l) = Tj&ftA; {wf s) } L ) + 0(e) , (3.21) 



where T^.^li^ {w^ zs> }l) commutes with the projection operator P^'l as follows [T5] . 

p(2s) T (l,l+) n r (2s)-, n_ p (2s) T (l,l+),,.r (2s)-, . p (2s) 
■ / 12...L J 0,12-£l A > i W j IL) ~ ^12 -L J 0,12-lI A ' \ W j h)^12-L ■ 

We show <K2H\f making use of fl2~22|) . 



,(28) 



(2s) 



(3.22) 
□ 
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3.3 Expectation value of a local operator through the limit: e — > 

In the massless regime, we define the expectation value of product of operators rifcLi E % ^^ k <y2s+ ' ) 

— ^(2s+) 

with respect to an eigenstate |{A Q } M ) by 

/TT p i fc.J'fc( 2s +)\ ( J\ \( 2 «+)\ _ \I A a/Af I ilfc=lj^fc ll A Q/Af / /o O o\ 

(11% / (^ivu J- — (2a +) ^^(2s+) • y 6 - 26 ) 

k=1 ({K} M \{K} M ) 

We evaluate the expectation value of a given spin-s local operator for a Bethe-ansatz eigen- 
state |{A a }^), as follows. We first assume that the Bethe roots {A q (c)}a/ are continuous with 
respect to small parameter e. We express the spin-s local operator in terms of spin- 1/2 global 
operators such as formula (13. 18j) with generic inhomogeneous parameters w^ s,t \ Applying 
(13.191) and (13. 19j) we remove the projection operators out of the product of global operators. 
We next calculate the scalar product for the Bethe state ({A^e)}^'^) which has the same 
inhomogeneous parameters w^ 2s ' e \ making use of the formulas of the spin-1/2 case. Then we 
take the limit of sending e to 0, and obtain the expectation value of the spin-s local operator. 

For an illustration, let us consider the expectation value of _g™' n ( 2s +)_ First, applying 

projection operator P^ s l L to |{A a } M ) = n^Li ^ 2s+ \X a )\0) we show 

-(2,+), 



^ ( . 2 S|{AJm > = A (2 2n fi(2S+;e) (A«(e))|0) + O(e) 

a=l 

M 

= e-^- A ^) P?ixn...L n 5(2s;e) (Aa(e))|0) + O(e). (3.24) 



a=l 

Second, making use of the relation (0| = {0\P^ S } L , we show 
1 - s + 1 



({AJm^I = (0| nC (2s+;e) (A Q (6))Pffi + 0(6) 

a=l 
M 

- (0| n^^Aa^xi-x^^^ + OCe). (3-25) 



a=l 



Making use of (13.181) we have 

^-^(2s+) ~ nn (2s+) -~^( 2s +) 

({A Q } M \E" n s |{A Q } M ) 



2s 
n 



M n 2s 



(oin^ (2s+; ^A a ( e ))ps^n^ (2s+;e) (^ ;e) ) n 

a=l k=l k=n+l 

M 

11 (A( 2s+ ->^ + D^ 2a+ ^)(w^) P^l- Y[B (2s+) (\ a )\0) + 0(e). (3.26) 



Here we have n 2 =L S {A( 2s+ ' ^ + D^ 2s+ ' , ^){wf s,e " > ) = I® L for generic e. We apply projection 
operators p( 2s ) to _p( 2s ) from the left, which are underlined in (13.261) . and make use of (I2.22p . 
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We then move the projection operators p( 2s ) in the leftward direction, making use of (I3.19p . 
Thus, the right-hand side of (13.261) is now given by the following: 



2s 
n 



M 

(2«+;e)/,.,(2»;eh 



(oin^ (2s+;e) (A«( e )) n^ (2s+;e) ^f ;e) ) n Ai2 

a=l k=l k=n+l 

M 

Yl (A {2s+ ' e) + D {2s+ ' €) )(wf s ' e) ) f]5 (2s+;e) (A /3 (e))|0) +0(e). (3.27) 



After applying the gauge transformation Xi~-l i nverse to (12. 7p [13] . we obtain 

-—^ — (2s +) ~ ,„ ,» - — - — (2s +) 

„ \ M n 2s 

iim(oin^ ;e) (A Q (e)) n^ (2s;e) (^ 2s;e) ) n A™{w™) 

U J € q=1 k=l k=n+l 

2sN 3 M 

x Y[ {A^ + D^)(wf S]e) ) Y[B^ e \X^e))\0) . (3.28) 

j=2s+l (3=1 

Here y4^ 2s;e ^ and _D( 2s;e ) denote matrix elements A^ 2sp '^ and D^ 2sp]e ^ of the monodromy matrix 
with principal grading, respectively. In the last line of (I3.27p . we have evaluated the eigenvalue 
of transfer matrix A^^lw^'^) + D^ 2s ' ,e \wf s ' ) on the eigenstate |{A^(e)}^- s '^) as follows. 

2sN s 



J] {A^iwf^) + D^\wf^)) \{X,(e)}^) 

j=2s+l 

2sN s M \ 

n n ri ( A ^)-^ 2s;e) ) i{A^)}& ;e) )- (3.29) 

\j=2s+l a=l / 

Before sending e to 0, we expand the products of C operators multiplied by operators A 
and D by the commutation relations between C and A as well as C and D, respectively. We 
then evaluate the scalar product of B and C operators with inhomogeneous parameters wf s ' £ \ 
Finally, we derive the expectation value in the limit of sending e to 0. 

Sending e to 0, we calculate the expectation value of A^^X) + D^(X) at A = wf s \ For 
instance, we calculate A^' e \w^ s '' e) ) + D {2s ^\ w^ s ' e ^) on the vacuum |0) as follows. 



lim(0| [A^\wf s ' e) )+D^\wf s ;e) )J |0) 

iim(o| (A(wf' e) ; {wf s ' e) } L ) + D^(wt' £) ; {w?*>} L j) |0) 

/ tN. \ 

Jim 1 + ]Jb(wt ;e) - wf s;e) ) (0|0) 

(1 + 0)(0|0). (3.30) 
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If we put A = after sending e to 0, the result is different from (I3.30P as follows. 



lim (0| (>*)(A; {wf s) } L ) + D ( - 2s \X; {wf s } } L )))|0) = 1 + \[h(wf s ) - Q (0|0) . 

X ^ W 2 \ p=l J 

(3.31) 

4 Derivation of matrix S 

4.1 The ground-state solution of 2s-strings 

We shall introduce ^-strings for an integer £. Let us shift rapidities Xj by sr] such as Xj = Xj + srj. 
Then, the Bethe ansatz equations (13.91) are given by 

ff Sinh( ^ Z*p± ^ = ff ^i-h + v) t fatj = li2l ...,n. (4.1) 
smh(A i - C P - sr]) pJi^ smh(Aj -Xp-rj) 

We define an ^-string by the following set of rapidities. 

Ai a) =/ia + (^+l-2«)| + e^ for a = l,2,...,£. (4.2) 

We call fi a the center of the ^-string and eh* 1 string deviations. We assume that eh* 1 are very 
small for large N s : 

lim e M = 0. (4.3) 

If they are zero, then we call the set of rapidities of (14.21) a complete ^-string. The string 
center /i a corresponds to the central position among the £ complex numbers: A^, X^\ • • • , \^P- 
Furthermore we assume that ]i a are real. If inhomogeneous parameters, £ p , are small enough, 
then the Bethe ansatz equations should have an ^-strings as a solution. 

In terms of rapidities Xj which are not shifted, an ^-string is expressed in the following form: 

A^ = ]i a - (a - 1/2)77 + 4°° for a = l,2,...,£. (4.4) 

We denote X^ also by A( a)a ). 

Let us now introduce the conjecture that the ground state of the spin-s case \ipg) is given 
by N s /2 sets of 2s-strings: 

A^ = fx a — (a — l/2)?7 + 4 a) , for a = l,2,...,N s /2 and a = 1,2,..., 2s. (4.5) 

In terms of Aa s in the massless regime, for w = + and p, we have 

N s /2 2s 
a=l o=l 

Hereafter we set M = 2sN s /2 = sN s . 
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According to analytic and numerical studies [101 SH E2] , we may assume the following 
properties of string deviations e„ s. When N s is very large, the deviations are given by 

(4.7) 

where i denotes \/—T, and 5^ are real. Moreover, 5a — 5 a a > for a = 1, 2, . . . , 2s — 1, 
and \5^\ > \5^ +1 ^\ for a < s, while \5 a \ < |<5a* +1 ^| for a > s. 

In the thermodynamic limit: N s — > oo, the Bethe ansatz equations for the ground state 
of the higher-spin XXZ chain become the integral equation for the string centers, as shown in 
Appendix F [53J. The density of string centers, p tot (/i), is given by 

Ptot(/i) = 77 y^T^T tj—( ~ t \ it\ (4-8) 

Ns ~[ 2 C cosh(7r(/i - Q/Q 

Thus, the sum over all the Bethe roots of the ground state is evaluated by integrals in the 
thermodynamic limit, N s — > oo, as follows. 

M , 2s N s /2 



s A=l s a=l a=l 

= E / ^ " ( Q " ^ + ^ W d ^ + • (4-9) 

a=l 

For the homogeneous chain where £ p = for p = 1, 2, . . . , N s , we denote the density of string 
centers by p(A). 

* A > = ycJfrVfl - (4 ' 10) 

Let us introduce useful notation of the suffix of rapidities. For rapidities = \(a,a) we 
define integers A by A = 2s(a — 1) + a for a = 1, 2, ... , iV s /2 and for a = 1, 2, ... , 2s. We thus 
denote A( a , a ) also by A^ for A = 1, 2, ... , siV s , and put A( a)Qi ) in increasing order with respect 
to A = 2s(a - 1) + a such as A^i) = Ai, A( lj2 ) = A 2 , • • • , X(n s /2,2s) = X sNs . 

In the ground state rapidities A^ for A = 1, 2, ... , M, are now expressed by 

A 2s( a-i)+a = Ma - (a - 1/2)77 + 4 Q) for a = 1, 2, ... , iV s /2 and a = 1, 2, ... , 2s. (4.11) 

For a given real number x, let us denote by [x] the greatest integer less than or equal to x. 
When A = 2s(a — 1) + a with 1 < a < 2s, integer a is given by a = [(A — l)/2s] + 1, and 
integer a is given by a = A — 2s[(A — l)/2s]. 

4.2 Derivation of the spin-s EFP for a finite chain 

We define the emptiness formation probability (EFP) for the spin-s case by 

{ 2 S+)(m) (4 2S+) |^' 2S(2S+) ---^' 2S(2S+ ^f +) ) 

(4 2s+) i4 2s+) ) 



19 



We shall denote T^ 2s+ \m) by r^ 2s ^(m). 

Let us assume that Bethe roots {A q (c)}m with inhomogeneous parameters wf s,e ^ (j = 
1,2, ... ,L; L = 2sN s ) become the ground-state solution of the spin-s XXZ spin chain, {X a }M, 
in the limit of sending e to 0. We denote the Bethe vector with Bethe roots {X a (e)}M by 

M M 

l^ 2S+;£) ) = n 5(2S;e) ( A ^ e ))l°) = e " E " ilAQ(e) ^...L ■ l[B^\X a (e))\0) 

= e-^ x ^ Xl2 ... L \^). (4.13) 

Here we recall the transformation inverse to (12. 7p . We now calculate the norm of the spin-s 
ground state from that of the spin-1/2 case through the limit of sending e to as follows. 

(^f +) l^f +) ) = lim(< 2 -)|^ 2 -)) 

M M 

= Jim(0|nC' (2s;e) (A fc )n s(2s;6) (Ai)|0) 
e ^ k=i j=i 

M 

= limsmh M ry J] ^(A^e), A*(e)) ■ det<f>«' ({A fc (e)} A/ ; {wf 8 

j,k=l;j^k 
M 

= sinh A S J] b-\\ s ,\ h ) ■ det$^'({A fc } M ;{£p}ivJ (4.14) 

3,k=l;j=£k 

where matrix elements of the spin-s Gaudin matrix for j, k = 1, 2, . . . , M, are given by 

*g? '({*}*;{&}) 

d , (a^(Xj) rr sinh(A t -A J +r / )' 



«9A fc ta \^ 2s )(A,) si nh(A 4 - A, - n) 
( smh.{2srj) ^— > sinh2^ 



-f sinh(Aj - £ p ) sinh(Aj - £ p + 2sr/) ' sinh(Aj - Ac + r]) sinh(Aj — Ac — rj) 

^ sinh2?7 ^ 

sinh(Aj — Afc + 77) sinh(Aj — A^ — rf) 

By applying formula (I3.18P with n = 2s, the numerator of (14.121) is given by 

m 

^s + ;,)^2s,2s(2s+) ^ 2s (2s +) ^(2s +; «)) = j^( 2s; e) | "Q ^ )2s ( 2s ) ^( 2s;e)) 

fc=l 

m /2s(i-l) 2s 



iim(^ 2s;e) i4 2 i n n + ^ (2s;e) ) (^ 2s;e) ) • n^^av*) 
•n^ (2s;e) +^ (2s;e) )(^ 2s;e) ) 4 2 ii^ s) ) 

a=l J 

m M 

nn & -(A«,0) Km^l ^ (2s;e) rf S;£) ) • ••/^-(- 2 ;;; : ) |^ 2s;e) ) • (4.16) 
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Let us set AM+j(e) = w f S '^ f° r 3 — 1? 2, . . . , 2sm. Applying formula flG.lj) to (14.161) we have 

M 2sm M 

(0|n^ (2s;e) (A a (e))n^ (2s;e) (AM +J (e)) J] B^A^e)) |0> 

o=l 3=1 0=1 

MM M 

= Yl Yl " C Cl ... C2sm (Ai(e),--- ,A M + 2s m(e);{^f s;e) }L) 

ci=l c 2 =l;c 2 ^ci C2 sm =l;C2 sm ^ci,...,c 2sm _i 
M+2sm M 

x(0| J] C( 2s ^)(A fc (6)) n^ (2s;e) M 2s;e) )|0), (4-17) 

k=l;k^Cl,...,C2sm Q=1 

where 

G Cl ... C2sm (Ai, • ■ ■ , A M +2 S m; Wli) = ]M a(A Cj ; {wj-J-l) — — jg+j — — — — 

j=i \ lli=i;t^ci,..., Cj . smn(A Cj . - \ t ) 

(4.18) 

We remark that from (I4.18P the set of integers ci, . . . , C2 Sm of the most dominant terms in 
(14.171) are given by m sets of 2s-strings. If they are not, the numerator of (I4.18P and hence the 
right-hand-side of (I4.17P becomes smaller at least by the order of 1/N S in the large N s limit. 
However, each of the most dominant terms diverges with respect to N s in the large- N s limit, 
and they should cancel each other so that the final result becomes finite. We therefore calculate 
all possible contributions with respect to the set of integers, ci, C2, . . . , C2 Sm . 

Let us take a sequence of distinct integers Cj satisfying 1 < c 3 < M for j — 1, 2, . . . , 2sm. 
We denote it by (cj)2 Sm , i.e. (cj)2 sm = (ci, C2, . . . , C2 Sm ). Let us denote by Em the set of 
integers, 1,2,..., M: S M = {1,2,..., M}. We then consider the complementary set of integers 
S M \ {ci, . . . , c 2sm }, and put the elements in increasing order such as Z\ < z 2 < ■ ■ ■ < z M _2 Sm - 
We then extend the sequence z n of M — 2sm integers into that of M integers by setting 
Zj + M-2sm = Cj for j = 1, 2, ... , 2sm. We shall denote z n also by z{n) for n = 1, 2, . . . , M. 

In terms of sequence (z n ) m we express the scalar product in the last line of (14.171) as follows. 

M+2sm M 

<0| n C^\X k (e)) l[B^\X a (e))\0) 

fc=l;fc^ci,...,c 2sm a=l 

M—2sm 2sm M—2sm 2sm 

= (0| n C^\X z(k) (e))]JC^\wf^) B^\X z{l) (e))]jB^(X C] (e))\0). 

k=l j=l i=l j=l 

(4.19) 

We evaluate scalar product (14. 19j) . sending v 3 - to \ z (j)(e) for j < M — 2sm and to wfl'M+2sm 
for j > M — 2sm in the following matrix: 

# ((Az(fc)(e))M, { v z{l)i • • • ) V z (M-2sm), ^M-2sm+l, ■ ■ ■ , ^m)] {w r e) ) L ). (4.20) 
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Here we define the matrix elements H^ b s \{X a } n , {fij} n ] {£k}N s ) for a, b = 1, 2, ... , n, by 



- (^b Ji sinh(A * - * + - JL smh(Afc - * - 

(4.21) 

Let us denote M — 2sm by M . We write the composite of two sequences (a(i))M and 
(b(J)) N as (a(i)) M #(b(j))N- Explicitly we have 

K*)W#(&0'))jv = (a(l), . . . , a(M), 6(1), . . . , 6(N)) . (4.22) 

For j > M' = M — 2sm, we have 

,(2s;e)-. v . / NX f-ii \ .411,, .\ . f„,P^^h 



N s 



j-M 



' sinhr; 



II sinh(A a (e) - w { *$ + rj) 



a=l 



' u ' » sinh(A 2(i) (e) - w^'fy smh(A z(i) (e) - w^] + rj) 



^-"" { ' >l) fisinh(A ((£ )- ro f-; + ,) 

sinh?7 



(2s; e) v. . w, , \ (2s; e) 



sinh(A zW (e) - ur.^J sinh(A z(l) (e) - + 77) 



(4.23) 



The second term of (14.231) for matrix element vanishes since we have d(w {w^ }l) = 
0. Here we remark that if we directly evaluate matrix if ( 2s ) at e = 0, the second term of (14.231) 
for matrix element for j 7^ 2s(n — 1) + 1 + M' with n = 1,2, ... , m, does not vanish, 

although it is deleted by subtracting column j by column j — 1, as discussed for the XXX case 
in Ref. |13|. We thus have 



limdet# (1) ((A 2(fe) (e))A/, (A z( i)(e), . . . , X z{M -2sm)(t), w ( ^ s '' e \ . . . , wf s s ^ ] ); (wf S]e) ) 2sNs ) 

M-2sm M 2sm M 

= ( _ 1) M-2s TO "Q Y[smh(X k -X Zb -r ] )l[l[smh(X k -wf s) + V ) 

6=1 fe=i i=l fc=i 

x det¥ 2s y((X z(l )) M , (X z(i )) M/ #(wf s) ) 2sm , (e P K) (4.24) 

where element of ^ (2s) '((A 2 ( fc) ) M , (A 2 ( fe )) M '#(wf s) ) 2 s m ; {^ P }vJ for i = 1,2, ...,M, are 

given by 

(('W)' ■ ■ ■ J ^«(M-2sm), A C1 , . . . , A C2sm ), (A z (i), . . . , X z ( M -2sm)-,wf S \ . . . , «4ot»); (£p)iV 8 ) 

' 3^, W ((Aj0i#;{&}) for J<M-2,m, 

sinhr] . . , , (4.25) 

for j > M - 2sm . 



k sinh(A z(i) - tw]_'^) sinh(A 2(i) - u/_ s ^, + 77) 
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Therefore, for i,j = 1, 2, . . . , M — 2sm, we have 

($ (2s) '((A,( fc) )M; ^ (2s) '((A 2(fc) )M, {K{k)) M >#{™f S \sm] {£ P })) = kj • (4-26) 



In terms of sequence (cj)2 Sm , we express the dependence of matrix (<£>( 2s )') 1 v[/( 2s )' on the 
sequence of Bethe roots (\ z (i))M etc., briefly, as follows. 



(<S>W)-^y(( Cj ) 2sm , {Q) = (<&<*>')-i*<*>' ((A 2(J) ) M , (A, (i) V#( W f } ) 2 . m ; {&}) • (4-27) 
Recall M' — M — 2sm. Similarly, we define $( 2s ) '((c i ) 2sm , {£ P }) and * (2s) '((c 3 -) 2sm , {£ p }) by 

<f ((c fc ) 2sm ,{U) = ^ (2 ; ) '((A z(fe ))M,(A 2(fe )) A/ #( W f ) )2s m ;fe)7vJ, (4.28) 

for i,j = 1,2, ... , M. Here we remark again that sequence (z(i))m is determined by sequence 
(ci, . . . , c 2s m) by the definition that {z(l), z(2), z(M')} = {1,2, ... ,M} \ {ci, . . . , c 2sm }, 
and z(l) < - - - < z(M') while z(j + M') — Cj for j = 1,2,..., 2sm. 

From property (I4.26P we define a 2sm-by-2sm matrix < ' 2s;m - ) ((c :) ) 2sm ; {£ p }) by 

(2s;m) ((c,W; {^}) iiJfe = (($( 2s )') _1 ^ 2s )'(( Cj ) 2sm ; {Q)) for j, fc = 1, 2, ... , 2sm. 

V V 7 J j+M' ,k+M' 

(4.29) 

Making use of flG.2j) . we obtain the spin-s EFP for the finite-size chain as follows. 



r ^\ m ) = rr „:^m, 7v~ x 



ni< K ,< 2s smh m (r-j> " ni< fe <K m n?i 1 niisinh(e fc -6+(r-j» 

A/ M M 

X E E E ^ (2s HA Cl ,---,A C2sm ;{U)det(^ m )((c,) 2sm ;{U)) 

Ci=l C2 = l;C27^Ci C2sm = l;C2sm : ^Ci,...,C2 S m-l 



(4.30) 



2sm m 2s— 1 



where if (2s) (A Cl , . . . , A C2sm ; {£ p }) is given by 
iJ^(A Cl ,...,A C23m ;{U) 

= n sinh ( A -a +7?) x n n n - & + m 

m 2 s /' — 1 m \ 

x n n n ^ - & + n ^ - & ) ■ ^- 3i ) 

i=l rj=l \6=1 fc=«+l / 

4.3 Diagonal elements of the spin-s Gaudin matrix 

Let us define K n (X) for 77 = i( with < ( < ti by 

XJX) = — sinh(m/) 

27ri sinh(A — nr]/2) sinh(A + rvq/2) 
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Lemma 4.1. For < ( < n/2s, we have 

/oo 

K 2 (X- p + nr] + ie)p(p)dp (n = 1, 2, . . . , 2s - 1) , (4.33) 
-oo 

and 

poo 

K 1 (X-nr])= K 2 (X- p-nrj-ie)p(p)dp (n = 1, 2, . . . , 2s - 1). (4.34) 



Here we recall rj = i(. 

Proof. We first consider the case of positive n. Let us recall the Lieb equation 

/oo 
K 2 {X-p)p{p)dp. (4.35) 
-oo 

Shifting variable A analytically to A + i( — ie in (I4.35P we have 

/oo 
K 2 (X — p + i( — ie)p(p)dp 
-oo 

Using 



(4.36) 



' 1 +2m5{p:-X) (4.37) 



sinh(A — p — ie) sinh(A — p + ie) 
we have 

poo POO 

K 2 (X- p + i(-ie)p(p)dp= / K 2 (X- p + i( + ie)p(p)dp + p(X). (4.38) 

J —oo 

Combining p(X + i() = —p(X) we obtain eq. (I4.33[) for n — 1. Making analytic continuation 
with respect to A we derive eq. (I4.33P for n — 2,3, ... ,2s — 1. Similarly, we can show (14.341) . □ 

Proposition 4.2. When < £ < ir/2s, matrix elements {A, A) of the spins Gaudin matrix 
with A = 2s (a — 1) + a are evaluated by 

' ^aa({^}m) = PtoM + 0(1/N S ) . (4.39) 



2m N, 

Relations ^4-39 ) are expressed in terms of integrals as follows. 



Ns 2s i*oo 

Ptot(Pa) = -Tf^2K 2s (p a -(a-l/2)f]-^ p +sr])-^2 / K 2 (p a -p c -{a~^)ri+e {an) )p tot (p c )dp c , 

s p=l 7=1 

(4.40) 

where e^ 1 " 1 = e^ — 6c . We recall that C corresponds to (0,7) with C = 2s(c — 1) + 7. 
Proof. Let us first show 

2s /»oo 

K 2s (p a -(a- 1/2)7] + S r])-J2 K ^a -Pc+(j- a) V + e {a ^)p(p c )dp c = p(p a ) (4.41) 

J -co 
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for a = 1,2 ... ,2s. Making use of the following relations 

K (X) = — ( cosh ( A ~^/ 2 ) _ cosh(A + m?/2) \ 
n{ ' 2m Vsinh(A -rvq/2) sinh(A + nrj/2) J { ' ' 



we have 

2s-l 



K 2s (\ + (2s - 1)^/2) = ^ i^(A + nr,) . (4.43) 



n=0 



We thus obtain (14.411) as follows. 



2s /»oo 

^2s(/ia ~(a- 1/2)7] +ST])-J2 / ^2(/ia ~ Pc + (7 ~ + ^"'^ )p(Pc)d/J c 

J= l 

2s / «oO \ 

£ (tfi(// a + (7 - a)r/) - J K 2 (fi a - (i c + (j - a)rj + e^)p(fx c )d^ c j 



= Kx^a) - I K 2 (p a - /J, c )p(^ c )dfjL c 
J — oo 

= P^a). (4.44) 

Here, in the second line of (14.441) . the summands for 7^a vanish due to lemma |4~T1 We then 
apply the Lieb equation (H~35) to show the last line of (T4T4l|) . We obtain (H~4"UI) from flOTj) . □ 

Corollary 4.3. Let us take a sequence of integers, C\, c 2 , . . . , c 2sm , which satisfy 1 < Cj < M for 
j = 1,2, . . . , 2sm, and determine a sequence (z(n))M by the conditions that {^(1), z(2), . . . , z(M')} 
{l,2,...,M}\{c 1 ,...,c 2sm } ! withz(l) < ■■■ < z(M') andz(j + M') = cj for j = 1,2, ...,2sm. 
Here we recall M' = M — 2sm. In the region: < £ < n/2s, diagonal elements of the 

spin-s Gaudin matrix $( 2s )' \(ck) 2sm ) are evaluated as 



1 

2m N, h0 



$7 {{c k ) 2sm )= p,ot{Pa) + 0{l/N s ), for 3 = 1,2,..., M. (4.45) 



Here integer a satisfies z(j) = 2s (a — 1) + a for an integer a with 1 < a < 2s. 



4.4 Integral equations 

We calculate matrix elements of (($( 2s )') _1 \J/( 2s )') ((cj) 2sm ) through the spin-s Gaudin matrix. 
For j,k = 1,2, ... , M, we have 

h i2sY (( Cj ) 2sm )) = Uw^'rvy) 

\ / j,k \ / j,k 

M 

= £ (^'(kW)) . ((* (2s) y^ (2s) '((9)2 Sm )) ■ (4.46) 
t=l *' 

We remark that matrix elements (A, B) of \I>( 2s )' \(ci) 2sm ) with A = j + M' and B = k + M are 
expressed in terms of i^i(A) as 

^ U+ M'(( c ^2 S m)/27rz = K 1 (A c .- U ;f ) + 7 7 /2) for j, k = 1,2, . . . ,2sm. (4.47) 
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Suppose that we have a sequence (z(n)) M for a given sequence (cj) 2sm satisfying 1 < < M 
for z = 1,2, ... , 2sm. Let us take a pair of integers j, k with 1 < j, k < M. We denote z(j) 
by A, and we introduce a and a by A = 2s(a — 1) + a with 1 < a < N s /2 and 1 < a < 2s. 
Applying proposition 14.21 and corollary 14.31 to (14.461) we have 



M 



i(cihs m )/2m ($(-)')-^(-)'(( 



t=i 



\Cl)2sm, 



t,k 



smh{2sr]) 



g{ (g i sinh(A A - Jsinh(A A - £ p + 2^) " g ^ " Ad) ) ^ 



*(*) 



+ if 2 (A A - A, (0 ) 1 (($ (2s) ')-^^'((Q) 2sm )) 



iV s Ptot(^)(($ (2s) ')- 1 ^ (2s) '((c i ) 2sm )). +^^ 2 (AA-A, {t) ) (($ (2s) yV 2 *)'(( Q ) 2sm )) +0(l/iV s ) 

(4- 



t=i 



Let us discuss the order of magnitude of the correction term in (I4.48p . It follows from 
( 14.451) that if the density of string centers p(p a ) is 0(1) in the large N s limit, then the diagonal 
element (A, A) of $( 2s )' is 0(N S ). We thus suggest that the matrix elements of ( - 2s;m ' ) ((Q) 2sm ) 
should be at most of the order of 1/N S , and hence the correction term in (I4.48P should be at 
most 0(1/N S ). 

Let us now define (Pa,b({£, p }) by the following relations for j, k — 1, 2, ... , M: 

¥>.<*),* = ^.PtotW ((^ (2s) ')" 1 * (2s) '((Q) 2s m)) , (4.49) 

where integer a is given by a = [(z(j) — l)/2s] + 1. In terms of function a(z) = [(z — l)/2s] + 1 
we have 

M 

Y,K 2 {\a ~ \ z{t) ) (($ (2s) ')-^ (2s) '((Q) 2sm ) X 



A? 



t.k 



Vz{t), k 



= ^^ 2 (A A -A z(t) ) , 

2s N s /2 

= ^2j^^2(ptot(Pc)y 1 K 2 (X A -X {c>l) )if2s(c-l)+y,k- (4.50) 
7=1 s c=l 

Here we have replaced the sum over t by the sum over c and 7 where ^(t) = C = 2s(c — 1) + 7 
with 1 < 7 < 2s. Expressing z(j) and k by A and -B, respectively, we have the following 
equations: 

2s N s /2 

^M{U)+EFE^W) _1 ^( A ^^))^ B (fe}) = — ^(( Q ) 2sm )+o(i/iv s ). 

7 =1 s c =l 

(4.51) 
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Let us introduce b and (3 by k = 2s(b - 1) + (3 + M' with 1 < (3 < 2s and 1 < b < N s /2. In 
terms of string center /i a we express (or approximate) <p 2 (j),£:({£p}) by a continuous function of 
/i a and as follows. 

%(*),*({&}) = + C(l/iV s ) . (4.52) 

By taking the large- iV, limit, the discrete equations (I4.5ip are now expressed as follows. 

<P?(Pa, 6) + Yl / - + (7 - + ^c)<pf (p c , 6)^c = ^i(A c . , - ™f ) + v/2) . 

7=1 •'-oo 

(4.53) 

Here €ac = — We recall that for j > M' we have set z(j) = c-_ M > . 

Lemma 4.4. In the region < ( < tt/2s, a solution to the integral equations l[4-53\ ) for integers 
A = 2s(a - 1) + a (i. e. (a, a) ) and B = 2s(b - 1) + (3 + M' with 1 < a, (3 < 2s and 1 < b < m 
is given by 

<PA, B = (/ia, 6) = P(Pa - &) So, p. (4.54) 



Proof, (i) In (a, a) case, i.e. when integers A = 2s(a — 1) + a and B = 2s(b — 1) + a 
correspond to indices (a, a) and (6, a), respectively, assuming that (pj(/i c , £ b ) = for 7 7^ a, we 
reduce integral equations (I4.53P to the Lieb equation (I4.35p . Therefore, we have ip^a\fi a , £&) = 
p(p a -&)■ 

(ii) In (a,/?) case, i.e. when A = (a, a) and S = (b,(3) with (3 ^ a, assuming <Pt^(p c , £&) = 
for 7 7^ (3, we have from (I4.53P 

# 2 (p a - Pc + (P - + e AB ) V f\fi c , £ h )dfjL c = — ^% . (4.55) 

For a < (3, we have e^s = it. Shifting \i a analytically such — {(3 — a — 1)77, we have 



,09)/.. ^W.. 1 sinhr ? 



27ri sinh(/i a +i]-^ b -rj/2) sinh(/i a + 77 - £ 6 + tj/2) ' 

(4.56) 

Making use of (I4.37P we reduce it essentially to the Lieb equation. We thus obtain (p^(ji a , £ 6 ) = 
p(p a — £&). For a > (3, we have = — ^e, and show it similarly, shifting \i a analytically as 
Pa ^ Pa - (a - (3 + IK □ 

Proposition 4.5. Lei its ta£;e a set of integers, ci, . . . ,c 2sm , satisfying < Cj < M for j = 
1,2, . . . ,2sm. Suppose that the number of Cj which satisfy Cj — 2s[(cj — l)/2s] = a is given 
by m for each integer a satisfying 1 < a < 2s. Then, when < ( < tt/2s, the solution 
to integral equations \4-53\j for A = Cj with j = 1,2, .. . ,2sm and for B = B' + M' where 
B = 1, 2, . . . , 2sm, is given by 

vfK-,6) = PK- -Zb)6 a ,p, (4.57) 
where aj = [(cj - 1) /2s] + I, a = cj - 2s[{c j - I) /2s] and B' = 2s(b - 1) + (3 with \<(3<2s. 
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Proof. It follows from lemma l4~4l that <f C] ,B of (14.571) gives a solution to the integral equations. 
Taking the Fourier transform of (14.531) . we show in §5.5 that the set of integral equations (14.531) 
for A = Cj for j — 1,2, ... , 2sm and B — 1,2, ... , 2sm has a unique solution. Thus we obtain 
the unique solution (14.571) . □ 

Let us recall the assumption that function <pT\p> a , £&) is continuous with respect to fi a and 
Then, for any given set of integers, c%, . . . , C2 Sm satisfying < Cj < M for j — 1, 2, . . . , 2sm, 

we may approximate the matrix elements of ($( 2s )') _1 \[/( 2s )' as follows. For integers j and k 

with 1 < j,k < 2sm, we define bk and (3k as follows. 

a j = l( c j ~ 1 )/ 2s ] + l > a j = c j ~ 2s l( c j ~ l )l 2s \ , 

b k = [(k - l)/2s] + 1 , (3 k = k - 2s[{k - l)/2s] . (4.58) 

Then, we have 

((*(*)')-!*(*)' (( Ci ) 2sm )) =± Pi ^~^ k) S aj , Pk + 0{l/Nl) . (4.59) 

Here we recall M' = M — 2sm. 

For a given 2s-string, \{ a ,a)i with a = 1, 2, . . . , 2s, we define X^ aa ^ by the 'regular part' of 

A( a , a ): 

*'(«,«) = Ha ~ (a - 1/2)77- (4.60) 
Let us introduce a 2sm-by-2sm matrix S by 

Si,*(ci, • • • , c 2s m; (Qn s ) = p(^' Cj ~ wf s) + v/ 2 ) $aj,i3 k for j,k — 1, 2, ... , 2sm. (4.61) 
Here a,, aj, fe^ and (3k are given by (I4.58p . Then, we obtain 

^r\(c k ) 2sm , {&}) = 1 — J-^ % fc ((c fc ) 2sm ; (^)at s ) + 0(l/iV s 2 ) . (4.62) 



and we have 

2sm 



/ 1 1 \ 

d*(^((^; { «))=n(^^)-(d^(( 



c,- W; (£ p )atJ + 0(l/iV s ) . (4.63) 
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4.5 Fourier transform in the cases of spin-1 and general spin-s 

The integral equations (14.531) for the spin-1 case are given by 



<p { i ) (ji,Z P ) + I K 2 (fi-X) t p\ 1, (X,QdX+ I K 2 (ji- X + Tj + e^^iKQdX 

1 sinh(?7) 



,(1.2) VJI), 



(2) 



2m sinh(/i — £ p + |) sinh(/i — £ p — | ' 

/OO /*QO 
K 2 (/i - A)^S 2) (A,e p )^A + / K 2 (/i - A + + e^^A^dA 
-OO «/ — OO 



sinh(?7) 



2ni sinh(/x - £ p + |) sinh(/x - £ p + 



3r/ ■< 



^ 2 1] (M P ) + / ^ 2 (/x-A-r7W 2 ' 1) M 1) (A,£ p )dA + / K 2 {yL-XM\\QdX 

sinh(?7) 



2vri sinh( / u - £ P - f ) sinh(/x - £ p - \ 



/OO /'OO 
K2(/i _ A _ r/ + e (2,i) )v9 (2) (A ^ p)rfA+ / K% {p.-X)<pf{\QdX 
-OO J — OO 

1 sinh(^) 



27ri sinh(/x - £ p — |) sinh(/i - £ p + |) 



(4.64) 

We solve integral equations (14.531) via the Fourier transform. Let us express the Fourier 
transform of function (pffi (/x, £ ) by 



e^pj? (A*, for a, /3 = 1, 2, . . . , 2s . 



(4.65) 



We denote by K n {u) the Fourier transform of kernel K n (X). We define matrix by 



A* 



(2s) 



a/3 



ftJ?(u,S) for a, = 1, 2, 



2s. 



(4.66) 



We introduce a 2s-by-2s matrix M.% s ^ ■ We define matrix element (J, k) for j, k = 1, 2, . . . , 2sm, 
by 

r r°° 

1+ / e^ w K 2 (^)d^ for j = fc , 



3,k 



e ivu K 2 {[i + (k - j> + iO) dfi for j < k . 



(4.67) 



e VUJ K 2 (fi - (j - k)rj - iO) dfi for j > k. 
When < C < 7r/2s, we calculate the matrix elements of as follows. 



(2s) 
K 2 



5^(1 + K 2 {u)) + (1 - ^e^'^ (£ 2 (u;) - e^C?-*)^ f or j A; = 1,2,..., 2s. 

(4.68) 
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Here we define sgn(j — k) by the following: sgn(j — k) = — 1 for j — k < 0, and sgn(j — k) = +1 
for j — k > 0. Here, K 2 (lo), is given by 



,7T 



Ko(uj) 



e iuj ^K 2 (fi)dfi 



Similarly, we define a 2s-by-2s matrix by 



sinh(- - ()uj 
smh( — ) 



(4.69) 



M 



(2s) 



e^Kxifi - 6 + (Jfe - for j, fc = 1, 2, . . . , 2s. (4.70) 

When < Q < 7r/2s, we can show 

( M x?)j k = ^ {^lM + (1 - kk) e (k - j)C " {KM - e s ^- k ^ 2 ) } (4.71) 
for j,k = 1,2, ... , 2s. Here is given by 



Kt{oo) = / e^K^dfi 



sinh (|~^ 
smh(— ) 



(4.72) 



(4.73) 



Taking the Fourier transform of integral equations (14.531) we have the following matrix 
equation. 

m%>m™ = m%>. 

For the spin-1 case, from (14.641) we have 

1+K 2 (u) e<"K 2 (u)-l \ ( (p?(u) cpf\u) 

K x (u) e&Kxiw) - e^ 2 

e-^KM-e-W 2 K^uj) 

It is easy to show that matrix M^ 2 ' (tp) is given by the following: 

2cosh(Cu;/2) 




(4.74) 







\ 



(4.75) 



2cosh(Cw/2) / 



(2} 

We calculate the determinant of M.^ (the spin-1 case) as follows 



det 



l + K 2 (u) e< u K 2 (w)-l 
e-< u K 2 (w)-l 1 + K 2 (lu) 



1 + K 2 (lu) e< u K 2 (u)-l 
-(1 + e-^) 1 + e - Cw 

K 2 {u){l + e< u ) e^K 2 (cu)-l 
1 + 

K 2 (uj)(1 + e~^)(l + e^) . 



det 



det 



(4.76) 
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Here, we first subtract the 2nd row by the 1st row multiplied by e We next add the 2nd 
column to the 1st column. Finally, the determinant is factorized and we have the result. 

(2s) 

By the same method we can calculate the determinant of ^ or ^ ne s pi n " s case. We 

have 

det-Mg) = (2 cosh(C^/2)) 2s K 2s (uj) . (4.77) 

The determinant is nonzero generically, and hence the solution to matrix equation (I4.73P is 
unique. Therefore, we obtain the solution of integral equation (14.531) . 



5 The EFP of the spin-s XXZ spin chain near AF point 

5.1 Multiple- integral representations of the spin-s EFP 

Let us derive multiple-integral representations for the emptiness formation probability of the 
spin-s XXZ spin chain. We shall take the large N s limit of the EFP (14.301) for a finite-size 
system, and we replace rapidity X Cj with complex variable Xj for j — 1,2, ... , 2sm, as follows. 
For a given rapidity of 2s-string, Xa = /i a — (ct — 1/2)77 + e^, we define its regular part A^ by 
X' A = fi a — (a — l/2)r/. m the large- N s limit, we first replace A Cfe = X' Ck + e Cfc by X' k + e Cfc where X' k 
are complex integral variables corresponding to complete strings such as X' k = fi^ ~ (P ~ 1/2)77 
for some integer (3 with 1 < (3 < 2s where rj = i( with < £ < n and /!& is real. We express A^, 
and e Cfc simply by X k and e k , respectively, and then we obtain multiple-integral representations. 

Applying (I4.63P we derive the emptiness formation probability for arbitrary spin-s in the 
thermodynamic limit N s — > 00, as follows. 

r£ s \m- fe}) =n i < J<r < 2s (sinh(r- J »- X Ui<k<i<m Uf=i YllU sinh(& - & + (r - fjrjj 

2sm I 2s ^.oo+(-fc+i)r; \ 

x[] E dXl \ ^ (2s) (Ai, ■ ■ ■ , A 2sm ) det5'(A 1 , . . . , A 2sm ) 



=1 \fc=l J -oa+(-k+±)Ti 

where # (2s) ((-^)2sm) is given by 
# (2s) ((AzW ' 



(5.1) 



rii<i<fc<2 Sm sinh(A fc -X l + r ] + e k:l ) 

2sm m 2s— I 

x H ] [ I I sinh(A ; - £ fe + (2s - p) V ) 



(5.2) 



1=1 k=i p=i 

m 2s /I— 1 m 

x Yl II I JJsinh(A 2 s(i-i)+n - 6 + 2s?]) JJ sinh(A 2s (j-i) +n - 6 

1=1 n=l \k=l k=l+l 

and matrix elements of 5(Ai, . . . , A 2sm ) are given by 

Q J P( A i - 6 + {k ~ if ^ ~ H = (I - k )v o\ 

(J otherwise. 
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Here fij denotes the center of the 2s-string in which Xj is the kth rapidity. Explicitly, we have 
Xj = fij — {k — l/2)?7. m the denominator, we have set e^i associated with A^ and A; as follows. 

e _ I ie for lm(X k - A/) > , ^ 
I —ie for Xm(Afc — A;) < 0. 

In the homogeneous case we have e p = for p = 1,2, ...,N S . We have thus defined 
inhomogeneous parameters £ p . We recall that in the homogeneous case, the spin-s Hamiltonian 
is derived from the logarithmic derivative of the spin-s transfer matrix. 

Here we should remark that an expression of the matrix elements of S similar to (15. 3p has 
been given in eq. (6.14) of Ref. [13] for the correlation functions of the integrable spin-s XXX 
spin chain. 



5.2 Symmetric expression of the spin-s EFP 

We shall express the spin-s EFP (15.11) in a simpler way, making use of permutations of 2sm 
integers, 1, 2, ■ • ■ , 2sm, and the formula of the Cauchy determinant. 

Let us take a set of integers = a>2s(j-i)+k satisfying 1 < auj^ < N s /2 for j = 1, 2, . . . , m 
and k — 1,2, ... ,2s. Here we remark that indices d(j,k) correspond to the string centers ^(j,k) = 
[/,2s(j-i)+k- m order to reformulate the sum over integers, C\, . . . , C2 Sm , in eq. (I4.30p in terms of 
indices ^ us introduce c.\, . . . , c 2sm by 

c 2s (j-i)+k = 2s(a(j t k) - 1) + k , for j = 1,2, ...,m and k = 1, 2, . . . , 2s . (5.5) 

We also define (5{z) by (3{z) = z — 2s[(^ — l)/2s\. Then, dj are expressed as follows. 

cj = 2s{a j - 1) + (3{j) , for j = 1, 2, . . . , 2sm. (5.6) 

We decompose the sum over Cj into 2s sums over aj as follows. 

M 2s N s /2 

E 9{<H) = E E 4( 2s (% -!) + *)■ (5-7) 

Cj=l k = l CLj=l 

Let us consider such a function /(ci, C2, . . . , C2 S m) of sequence of integers (cj)2 S m that vanishes 
unless Cj's are distinct. We also assume that /(ci, C2, . . . , C2 Sm ) vanishes unless the number of 
Cj's satisfying /3(cj) = a is given by m for each integer a satisfying 1 < a < 2s. Here we 
recall that the two properties are in common with the summand of ()4.3Up . in particular, with 
det{^ 2s ' m \{ Cj ) 2sm ;{Q). Then, we have 

M M M N s /2 N s /2 N s /2 

EE - -- E f(°^--- > c 2*™) = r^yvTs E E ■■■ E E f^Ph---,c P{2sm) ) 

d=lc 2 =l c 2sm =l V ' a 1 =la 2 =l a 2 sm=i P&S2sm 

Ns/2 Ns/2 Ns/2 

= EE"' E E f(Cnl,--- ,CM2sm)) ■ (5.8) 

ai =la 2 = l a. 2s m=lTT£S 2sm /(Sm) 2s 
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Here an element 7r of <S 2sm /(<S m ) 2s gives a permutation of integers 1,2, ...,2sm, where irfs 
such that irj = k (mod 2s) are put in increasing order in the sequence (irl, ir2, . . . , 7r(2sm)) for 
k = 0,1,. . .,2s- 1. 

Reformulating the sum over Cj s in (14.301) in terms of a/s, in the large N s limit we have 
(2s)/ \ 1 1 

^ Ul<a< P <2s(^i0 ~ ^V) m X Ul<k<l< m II«=1 Iltl -& + (<*- 



X 



1/3= 

Is m /. oc j 

nn / d ^< k ) yi 7-^^ts{xp 1 ,...,xp(2s m ))H {2s \x P1 ,...,xp 



(2sm) , 



fc=i i=i " _CXJ Pe5 2sm 
1 

X 



ni< Q</3 < 2 >inh(/3 - a) V r Yli< k<l < m UlU njLi sinhfe - + (a - /fyj) 
x II / ^ 5Z detS , (A 7r i,...,A 7r(2sm ))iI (2s) (A 7r i,...,A 7r(2sm )), (5.9) 

3=1 °° 7r65 2sm /(5 m ) 2 = 



where symbols Aj denote the following 

A, = ^ - - ^ for j = 1, 2, ... , 2sm. (5.10) 
We calculate det5 applying the Cauchy determinant formula 

det f 1 \ = nr</ sinh (^ - 6) nr> fc smh(A Q - x b ) 
\sinh(A a -^)y nr=i nr=i sinh ( A a - 6) 

and we obtain the symmetric expression of the spin-s EFP as follows. 

T^(m-tt\) = - TT sinh 2 -(7r(q-a)/C) 

00 ^ ' UPIJ n x < a <,< 2s sinh m (/3 - a)n 11 itLr sinhfe - 6 + (r - J» 



x 



■2sm 2 /-■•-in ,. A \ 2- 



(2*0 



2sm 



11/ 1 1 1 1 sinh(7r(/i 2s(a _ 1)+7 - /i 2s (6-i)+ 7 )/C) 



VJ f=l " -00 / 7=1 i<6< a < m 
/ 2sm TT m tt2s— 1 



/ tt nr=i n/=i smh^ - & + ^ \ ^ 
x n nrioosh^.-^/o ^ (sgna) 

\ j=1 116=1 v vpj sb;/s; y <xeS 2sm /(S m )2s 

ni=i n 2 =i (nfci smh(A CT(2s( ,_ 1)+n) - ^ + 2s^) nr=«+i smh(A CT(2s(/ _i )+n) - 

v 

rii<i<fc<2sm sinh(A a(fc) - A ct(0 + 77 + e ff(fc)j(T(0 ) 

(5.12) 

Here (sgna) denotes the sign of permutation a G <S 2sm /(<S m ) 2s . 
5.3 The spin-s EFP for the homogeneous chain 

Sending to zero for 6 = 1, 2, . . . , m, we derive the spin-s EFP in the homogeneous limit. Here 
we remark that the spin-s Hamiltonian is derived from the logarithmic derivative of the spin-s 
transfer matrix t^ 2s ' 2s ^(A; {£b}jvj in the homogeneous case where = for b — 1, 2, . . . , N s . 
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Sending to zero for b = 1, 2, . . . , m, we have the following: 
lim Jim • • • lim (r^(m; {£ p } m )) 



sm(m-l) 



X 



X 



X 



ni< a</3 < 2 , sinh m {{p-a)ii) 

■2sm 2 "ism 2s 

-——J] / d^JJ J] sinh(7r(/i 2s(a _ 1)+/3 -^ 2s(6 _ 1)+/3 )/C) 
2 ^ n?= 1 1 sinh-(A J +^ 

n cosh-^/o J x £ ^) 

n^ Li rirjU (sinh'~ 1 (A g ( 2 5(i-i)+r,) + 2s??) sinh m ~'(A CT(2s(/ -i) +ri) )) 
rii</<fc<2 S m sinh(A ff ( fc ) - A ff (Q + r? + e CT (fc), CT (o) 



Here we recall definition (15. 10)) of Xj. 

Let us discuss that expression (15. 13)) gives the spin-s EFP for the homogeneous chain. First, 
we remark that T^ s \m; {£ p } m ) does not depend on £ p with p > m. Hence we may consider 
that inhomogeneous parameters £ p with p > m are all set to be zero, after computing the EFP: 

We now show that the order of the homogeneous limit: £ p — > and the thermodynamic 
limit N s —>■ oo can be reversed. We can show the following relation: 

m / \ ( m \ 

n*& K*V;{&u))=,ita n^ (^w&u) • (5.14) 

p=i v 7 \p=i / 

In fact, when N s is very large, it follows from (I4.63P that we have 

rt\m; {£ p } m ) = r^\m; {£ p } m ) + 0(1/N S ). (5.15) 

Furthermore, we can explicitly show that T^ s \m; {£ p } m ) is continuous with respect to £ p at 
£ p = for p = 1,2, ... ,m. We first reformulate the sum over Cj in (14.301) into the sum over 
a (j,k) by relation (15.81) . 

M M 

E d etS((c,) 2sm )tf^((A Cj ) 2sm ) 

Cl=l c 2sm = l 
m 2s f N B /2 \ 

= nnbE e dets((c PJ ) 2sm )^)((A £j ) 2sm ) 

j = l fe=l y ' a (j,fc) = 1 / fe&sm 
m 2s f N B /2 \ 

= nil La E ^(fo)*™) E (sgnP)i^((A % W). (5.16) 

j' = l fc = l \ ' 0(j iA )=ly Pe52sm 
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We then apply the Cauchy determinant formula to evaluate det S(dj) as follows. 

2s 

detS((c j ) 2sm ) = det5((A £; ) 2sm ) = JJ detS ,(a) ((A 2s ( a(j , a) -i)+a)m) 

0=1 

(n j<k sinh 40 - &vc) • ulu u j<k sinh {^(^s 



(2<)2 , m nr=i sinh ^ . _ & _ W2)/C 



(5.17) 



Making use of (I5.17P we show that such factors in the denominator of (I4.30p that vanish in the 
limit of sending £ p to zero are canceled by the factors in the numerator of (15.171) . We thus have 
shown that the EFP for the finite system, T^ s \m; {£ p } m ), is continuous with respect to £ p at 

e P = o. 

Therefore, expression (15.131) gives the spin-s EFP for the homogeneous chain. That is, we 
have the following equality: 

lim Km ... lim (ri 2s) K{U-)) = J™ r£ s) (m; {£ p = 0} Ns ) . (5.18) 

5.4 The spin-1 EFP with m = 1 

Let us calculate T^ 2s \m) for s = 1 and m = 1. From formula (15.11) we have 

1 / poo—'q/2 poo—Zri/2 \ / roo — 'q/2 poo— 3r//2 

r ( 2 )(l) = — — / dAi+ / dAi / dA 2 + / dA 2 

Sinn 7/ Y > /- o-r?/2 J-oo-3r)/2 J \J-co-r)/2 J-oo-3r?/2 

x^ 2 )(A 1 ,A 2 )det5(A 1 ,A 2 ) 

1 J_ r r d/i2 _ sinh( Ml - 6 + 77/2) sinh(// 2 - 6 - 77/2) 



i sin C 4C 2 7_oo 7_oo sinh^! -/i 2 + ie) cosh (tt^! -^)/C) cosh (vr(/i 2 -^)/C) 

1 (_i) p ,00 8^(^-^-77/2)8^2-^ + 77/2) 

a/ii / a/j 2 - 



i sin C 4C 2 7_oo i^oo sinh(//i - /i 2 - 2t?) cosh (7r(//i - f i)/C) cosh (7r(/i 2 - f i)/C) ' 

(5.19) 

Here we note that detS^Ai, A 2 ) = for Ai = \i\ —rj/2 and A 2 = /j 2 — 3r]/2, or for Ai = fi\ — 3i]/2 
and A 2 = /j 2 — 7]/2. Showing the following relations of integrals 

1 sinh(/j - 77/2) r , 1 sinh(/i + r//2) ,sinh(A - r//2) 

ct/i — — — . — — — I d/j — — — - r — 2ni- 



00 
00 



cosh (irfi/Q sinh(A — fi + it) J _ 00 cosh (irn/() sinh(A — fi — 77) cosh(7rA/£) 

dX 1 sinh(A + tt/2) = r , A 1 sinh(A-r//2) 

cosh (ttA/C) sinh(A — ji — 77) 7_ 00 cosh (nX/Q sinh(A — /U — 2rf) ' 



we thus have 



7T 1 / f°° cosh 2£x — cosh 77 



Evaluating the integral we obtain the spin-1 EFP with m = 1 as follows. 

r( 2 )(l) = C -; inC 2 rosC . (5.22) 
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Let us denote by (El' b ) the ground-state expectation value of operator E^' . For the spin-1 
case, we have E®' 0( ~ 2+ ^ + E^ 1 ^ 2 ^ + E 2,2 ^ 2 ^ = and hence we have 

(£0,0(2+)) + ( £M(2+) } + = i . (5.23) 

Due to the uniaxial symmetry we have (E^' ^ 2 ^} = (E 2 ' 2 ^ 2 ^} . Thus, we obtain 

( £M( 2+)) = cosC(sinC-CcosQ j (5 _ 24) 



lim S ZZ 2 T ( 5 " 25 ) 



In the XXX limit, we have 

C — sin C cos ( 1 
C-o 2( sin 2 C ~ 3 

The limiting value 1/3 coincides with the spin-1 XXX result obtained by Kitanine |13j . As 
pointed out in Ref. [13], (Ef 2 ) = (E\ l ) = (E°°) = 1/3 for the XXX case since it has the 
rotational symmetry. 

In the symmetric expression of the spin-1 EFP with m — 1, putting Ai = {i\ — r]/2 and 
^2 = ^2 — 377/2 in (15.131) . we directly obtain the following: 



T^(l) = ^—ip; I dfH I dfJL 



1 1 [°° f°° sinh(//i + 77/2) sinh(^ 2 - 77/2) 



zsinC4C 2 7_oo ^i-oo 2 cosh (vr/ii/C) cosh (7r/i 2 /C) 

x I — T7 r r? rl • (5-26) 

Vsmh(/i2 - /ii - le) smh(/ii - /i 2 + 2r]) J 

In the second line of (15.261) the first term corresponds to the first term of (I5.19p . while the 
second term to the second term of (15.191) with \i\ and /i 2 being exchanged. 



6 Spin-s XXZ correlation functions near AF point 

6.1 Finite-size correlation functions of the integrable spin-s XXZ 
spin chain 

We now calculate correlation functions other than EFP for the spin-s XXZ spin chain by the 
method of §3.3, making use of the formulas of Hermitian elementary matrices such as (13.181) . 

We define the correlation function of the spin-2s XXZ spin chain for a given product of 
(2s + 1) x (2s + 1) elementary matrices such as E^' h {2s+) ■ ■ ■ Em' jm {2s+) as follows. 

m 

F {2s+) ^ k ,j k }) = ^ +) i n^ fe,jfc(2s+) i^ 2s+) )/^f +) w 2s+) > . (6.1) 

k=l 

By the method of expressing spin-s local operators in terms of spin-1/2 global operators [15J, 
we express the mth product of (2s + 1) x (2s + 1) elementary matrices in terms of a 2smth 
product of 2 x 2 elementary matrices with entries {e^, e'j} as follows. 

m 2sm t 

JJ^,i*<*+) = C({i k ,j k };{e p e })P[^ L • • p£0.. (6 . 2 ) 

k=l k=l 
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We evaluate the spin-2s XXZ correlation function F^ 2s+ \{ik^ j^]) by 

2sm , 
3=1 

(6.3) 

We denote the right-hand side of (16.31) by F^ 2s+ \{ej, e'j}). 

Let us introduce some symbols. We denote by ck + the set of suffices j such that ej =0, and 
by a the set of suffices j such that €j — 1: 



0'; e i = o} , « = {j; e ',- = !} • 



(6.4) 



We denote by a + and a_ the number of elements of the set cx + and et~ , respectively. Due to 
charge conservation, we have 

a + + a_ = 2sm. (6.5) 

We denote by j min and j max the smallest element and the largest element of a~, respectively. We 
also denote by j' min and j' max the smallest element and the largest element of ck + , respectively. 

Let Cj (j G a~) and c'- (j G a + ) be integers such that 1 < Cj < M for j G a: and 
1 < c'j < M + j for j G ck + . We define sequence (fr^sm by 



(h,b 2 



' ^Jmin ? * * 



(6.6) 



Here sequence (c'., , . . . , c'., , Cj min , . . . , c-, max ) is given by the composite of sequence of c'-'s in 

3 max. -^min 

decreasing order with respect to suffix j, and sequence of c/s in increasing order with respect 
to suffix j. 

Let us introduce the following symbols: 











n 




n 


E 


ject- 









M 



M M+j min M+j n 



(6.7) 



r- — 1 ' 1 
L Jmax -i C ; =1 

7 . 
J mm 



Extending the derivation of the spin-s EFP we can rigorously derive the following expression 
of the spin-s XXZ correlation functions in the massless regime with < ( < n/2s. 









( ' M+j\ 






(!) 


n 


E 


detM^ 2sm \(h) 2sm ) 


jea- 











x(-iy 



Ujea (Ul=\ fi^ ~ wf s) + T]) nK+i ^(A c 



w 



(2s). 



Til<k<e<2sm f(\ ~ ^bh + V) 



n 



1 — w 



(2s) 



' — u» 



(2s). 



Ul<k<e<2sm ( P( wi k S) ~^? S) ) 



+ 0(l/iV s 



(6.* 
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Here y?(A) = sinhA. We have defined the 2sm x 2sm matrix M^ 2sm \ (pj)2sm) as follows. For 
£, k = 1, 2, . . . , 2sm, the matrix element of (£, k) is given by 

(M^\(b 3 ) 2sm )) p ={ ~ 6b r^ k foTb £ >M 

y U i> 2sm »^ ^ S m) , m ■ p(X k - wf s) + r 1 /2)/N s p tot {p a{bl) ) for bt < M 

where pj denote the centers of \j as follows. 

Aj = m - - \)r) 3 = 1, 2, • • • , 2sm. (6.10) 

We recall that a(j) and j3(j) have been defined in terms of the Gauss' symbol [•] by a(j) = 
[(j - l)/2s] + 1 and f3(j) = j - 2s[(j - l)/2s], respectively. 

We remark that under the limit of sending e to zero, the sum over variable cj is restricted 
up to M. 

6.2 Multiple-integral representations of the spin-s XXZ correlation 
function for an arbitrary product of elementary matrices 

Let us formulate matrix S for the correlation function of an arbitrary product of elementary 
matrices. We define the (j, k) element of matrix S = S (^(Xj) 2sm ; (w^^smj by 

Sj,k = P(A,-?4 2s) +r ) /2)5(a(X j ),P(k)), for j, k = 1, 2, . . . , 2sm . (6.11) 

Here 5(a, j3) denotes the Kronecker delta. We define a(Xj) by a(Xj) = 7 if Xj = Pj — (7 — 1/2)77 
or Xj = where f3(k) = 7 (1 < 7 < 2s). We remark that pj correspond to the centers of 
complete 2s-strings Xj. We also remark that the above definition of matrix 5* generalizes that 
of (15.31) since a(Xj) is now also defined also for Xj = w k 2s \ 

Let Tj be a small contour rotating counterclockwise around A = Wj. Since the detS has 
simple poles at A = with residue l/27rz, we therefore have 

/oo+ie poo—ie p 

detS((X k ) 2s m) dXi = detS((X k ) 2sm )dX 1 - ® detS((X k ) 2s m) dX x . (6.12) 

-oo+ie J —oo—ie ^Ti 

For sets ct~ and ct + we define Xj for j E a and A - for j E a. + , respectively, by the following 
sequence: 

(A' ,---,X, , X jmin ,Xj max ) — (Ai,..., A 2sm ) ■ (6.13) 

Jmax J rnin 

Thus, from the expression of the correlation function in terms of a finite sum (16.81) we obtain 
the multiple-integral representation as follows. 

F (2s+ \{e 3 ,e 3 }) 

oo+ie poo— i(2s— l)(^+ie \ / poo+ie poo— i(2s— l)f +ie \ 

+ •••+/ \dX 1 ■ ■ ■ / + •••+/ dX a+ 

-oo+it J —oo—i(2s—l)C+ie I \J —oo+ie J— oo— i(2s — l)(+ie J 

oo—ie poo—i(2s—l)<^—ie\ I poo—ie poo—i(2s—l)(—ie\ 

+ •••+/ dX a++1 ■ ■ ■{ +•••+/ ) dX m 

-oo—ie J —oo—i(2s—l)C,—ie I \J —oo—ie J — oo— i(2s— l)f — ie J 

xQ({e j , e^}; Ai, . . . , X 2sm ) det5(Ai, . . . , X 2sm ) ■ (6.14) 
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Here we have defined Q({ej, e' }; Ai, . . . , A 2sm ) in terms of small numbers of (15.41) by 



Q^e..}^,...^)) = (-!)«+ — V 



X u TIM (2sK (° 5 ) 



V - w k 



Thus, correlation functions (16. ip are expressed in the form of a single term of multiple integrals. 

Similarly as the symmetric spin-s EFP, we can show the symmetric expression of the 
multiple-integral representations of the spin-s correlation function as follows. 

F (*s+)(f e - e'\) = - TT sinh 2 -(vrfe-6)/C) 

U " ^ U^a^s sinh-(/3 - a) V J}^ ^i, UlU sinhfe - 6 + (r - J» 

°+ / /"oo+ie poo— i(2s— l)C,+ie \ 2sm / poo— it pco—i(2s — l)(—ie\ 

e n / +■■•+ / n / +■■•+ / 

o-G5 2sm /(5 m ) 2s i=l \ J -°°+K ./-oo-i(2s-l)C+ie/ j=a++l ^ -oo-i{2s-l)C-ie J 

/'S n^Li ri^i 1 s inh (^- - & + /fy) ' 

x (sgncr)g({e i ,e.};A ff i,...,A (7(2sm) )) II — = — — 

\f=T n&=i cosh(7r(/ij - 6)/C) 



•2sm 2 



X 



)/C) 



(2iC) 2sm 

v ^ 7=1 Kka<m 

It is straightforward to take the homogeneous limit: — > 0. Here we recall that (sgn a) denotes 
the sign of permutation a G <S 2 sm/ (<Sm) 2s ■ 

6.3 An example of the spin-1 correlation function 

Applying formula (13.181) to the spin-1 case with m = n = 1, we have 

2iV s 

£1,1(2+) = 2P^lD^+\ Wl )A^ + \w 2 )]J(A^ + D^)(w a ) P[ 2s) L . (6.16) 

0=3 

Therefore, we evaluate it sending e to zero, as follows. 
(^ 2+) l^' 1(2+) |^ 2+) )/(^ 2+) l^ 2+) ) 
= 2 lim(^ +;e) |i5 (2+;e) rf ;e) )A (2+;e) rf ;e) ) II( A(2+;£) + ^ (2+;e) )(^ 2;£) )l^ 2+;£) )/(^ 2+;£) l^ 2+;£) ) 

(poo+ie poo— if+ie\ / poo—ie poo— if— ie\ 

/ +/ ) dXi ( + )dX 2 Q{\ 1 ,X 2 )detS{\ 1 ,\ 2 ) (6.17) 

./— oo+ie J —oo—i^+ie/ \J—oo—ie J —oo—i^—it/ 

where Q(Ai, X 2 ) is given by 

<2 Ai, A 2 = (-1) tt 6.18 

y?(A 2 - Ai +77 + e 2> i)v9(r7) 
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and matrix S(\i, A 2 ) is given by 

s(x x)= ( P(*i-w[ 2) + v/m^iA) p(X 1 -w^+ V /2)6(X 1 ,2) \ 
1 2J \ p(X 2 -w^ ) +r ] /2)5(X 2 ,l) p(X 2 -w?+ V /2)5(X 2 ,2) J ' { ' ) 

We thus note note that the correlation function is now expressed in terms of a single product 
of the multiple-integral representation. 

Let us now evaluate the double integral (16.171) . explicitly. The integral over Ai is decomposed 
into the following: 

"oo+ie fao—i^+ie^ 



(roo+ie /•oo—zQ+ie\ 
/ +/ UA 1 Q(A 1 ,A 2 )detS(A 1 ,A 2 ) 
J —oo+ie J— oo— iC,-\-ie/ 

/ + / (iA 1 Q(A 1 ,A 2 )det5(A 1 ,A 2 ) 



- I dX 1 Q(X 1 , A 2 ) detS{\ u A 2 ) - / dX 1 Q(X 1 , A 2 ) detS(Ai, A 2 ) . (6.20) 
Thus, the integral (I6.17P is calculated as 

(^ 2+) i^ ii(2+ w 2+) )/f2*rw 2+ v 



OC 



. sinhfx — ri/2) sinhfx — 3ri/2) <,. . , , f°° sinhfx — ri/2) , , , 

-2m / i 'J—L- 5: l -^p 2 (x)dx + 2coshri / — -) !L-Lp( x )dx 

suahrj J_ OD smh(x + 77/2) 



f°° , x , / \ , f°° sinhfx - 77/2) . , , 

— / p(x)dx + ( — 1)2 cosh 77 / — — —p(x)dx 

J -00 K ' 1 ' / 7^ocSinh(x + 7 / /2) n ' 

cos C(sin C — C cos C) 

2Csin 2 C ' 

We have thus confirmed (I5.24p directly evaluating the integrals. 



(6.21) 



7 Concluding remarks 

In the paper we have explicitly shown the multiple-integral representation of the emptiness 
formation probability for the integrable spin-s XXZ spin chain in a region of the massless regime 
of 77 = iQ with < £ < -k /2s. We have also calculated the emptiness formation probability for 
the homogeneous case of the integrable spin-s XXZ spin chain. 

In the XXX limit where we send ( to zero, the expression of EFP for the spin-s XXZ case 
reduces to that of the spin-s XXX case. 

Moreover, we have presented a formula for the multiple-integral representation of the spin-s 
XXZ correlation function of an arbitrary product of elementary matrices in the massless regime 
where 77 = i( with < ( < tt/2s. We have also presented the symmetric expression of the 
multiple-integral representations of the spin-s XXZ correlation functions. 

Finally, we have introduced conjugate vectors ||2s, n) in order to formulate Hermitian ele- 
mentary matrices _g m > n ( 2s +) an d Hermitian projection operators P^' in the massless regime. 
We have also defined the massless fusion R- matrices R^' 2sw "> for w = + and p. 
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A Affine quantum group with homogeneous grading 

The affine quantum algebra U q (sl 2 ) is an associative algebra over C generated by Xf,Kf for 
i — 0, 1 with the following relations: 

K.K- 1 = K^Ki = l t KiXfK' 1 = q ±2 Xf , K % XfK^ = q ^Xf (i ? j) , 
[ X i~, X j] = kj q_ q -\ ' 

{XffXf - [3} q (Xt) 2 X±Xt + i3} q X±X±(Xt) 2 -X±(X±) 3 = (i ? j) . (A.l) 
Here the symbol [n] q denotes the g-integer of an integer n: 

n _ -n 

[n\ q = L_JL_ . (A.2) 

The algebra U q (sl 2 ) is also a Hopf algebra over C with comultiplication 

A(X+) = X+®l + Ki®X+, A(Xr)=Xr0K^ + l0Xr, 
A(Ki) = Ki^Ki, (A.3) 

and antipode: = Kr\S(Xi) = -K^Xf , S(Xr) = -XrR i: and counit: e(Xf) = 

and e(Ki) — 1 for i — 0, 1 . 

The algebra U q (sl 2 ) is given by the Hopf subalgebra of U q (sl 2 ) generated by Xf, Ki with 
either i — or i — 1. Hereafter we denote by X ± and K the generators of U q (sl 2 ). 

For a given complex number A we define a homomorphism of algebras (f\: U q (sl 2 ) — > U q (sl 2 ). 

ip x (X±) = e ±2X XT , Vx (Xf) = X ± , <p x (K ) = K- 1 , ip x {K x ) = K . (A.4) 

Map flA.41) is associated with homogeneous grading [8]. For a representation (n, V^) of U q (sl 2 ) 
we have a representation of U q (sl 2 ) by ir(ip x (a)) for a 6 U q {sl 2 ). We call it the spin-£/2 
evaluation representation with evaluation parameter A, and denote it by (7r A , V^(X)) or V^(A). 
We define opposite coproduct A op by 

A op (a) = ro A(a) for a G U q {sl 2 ) , (A.5) 

where r denotes the permutation operator: r(a ®b) = b ® a for a, 6 G U q (sl 2 ). 
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B Fusion projection operators being idempotent 

We give the derivation [23] of (-P^--^) 2 = P\2-i ma king use of the Yang-Baxter equations. 
Lemma B.l. Operators P12..4 defined by \2. 12\) have the following two expressions: 

Proof. Applying notation (12.21) to permutation operator we define permutation operators 
n^fc for integers j and k satisfying < j < k < L. The form of the left-hand side of (1B.1|) is 
expressed in terms of i?-matrices as follows (see also eq. (3.7) of [15J). 

[t/2] 

P\-t = II n i,^i+i • P-e-ie ' ' ' Rt,3-eRi,2-t ■ ( B - 2 ) 
3=1 

Making use of the Yang-Baxter equations (12.41) we reformulate (IB.2j) as follows. 

[1/2] 

P\-i = \\ n?,<?-i+l • Rt 2 Rl2,3 ' ' ' Pl,2-£ ( B - 3 ) 
3=1 

which gives the expression of the right-hand side of (IB.lj) . □ 

From fIB.ip we show that Pj^. 1 j +2 ...j + i is expressed as follows. 
p{t-i) R + ((P-l\ri\P il ~ l) -pC" 1 ) R+ ((£- l)n)p^ 

(B.4) 

Lemma B.2. Operator P^...^ projects operator Rf_ij(u) to 1 as follows. 

PiLRtiM = PiL- (B.5) 

Proof. Due to the spectral decomposition of the i?-matrix, we have P^ ie Rj~_ 1 e (u) = Pf\t 
Applying flBTTTl and flRlj) . we thus obtain (IB~5|) . □ 

Proposition B.3. Operators Pf^-i defined by h2.12\) are idempotent: {P^-.-ij = P\2-t ■ 
Proof. We show it from (1B.5|) by induction on I. Suppose that it is idempotent for i. We have 
(pd+i) V _ pit) B + (fri.fpM \2 .d+ (PripM 

- P (t} R + (fn) ■ P {e) H + (fn) ■ P {£) 

— r \ 2 -i ix ii^\y i, \) r i2---t JX it+i\ t 'l) r i2---e 

_ pit) d+ (f n \ pit) . pit) _ pit) p>+ if \ pit) (B fi \ 

□ 
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C Basis vectors of spin-^/2 representation of U q (sl2) 

In terms of the q- integer [n] q defined in ( 1A.2I) . we define the g-factorial [n] q \ for integers n by 

[n] q \ = [n} q [n-l] q ...[l] q . (C.l) 
For integers m and n satisfying m > n > we define the g-binomial coefficients as follows 

(C.2) 



m 
n 



TO 



[m - n] q \ [n] q \ ' 



We now define the basis vectors of the (£ + l)-dimensional irreducible representation of 
U q (sl 2 ), | \£, n) for n = 0, 1, . . . , t as follows. We define | \£, 0) by 



\\£,Q) = \0)i® |0> 2 <S> - - - <8> |0>^ . 



(C.3) 



Here \a)j for a = 0, 1 denote the basis vectors of the spin-1/2 representation defined on the jth 
position in the tensor product. We define \ \£,n) for n > 1 and evaluate them as follows 

1 



114 n) = {A^\X-)) n \\£,0) 



n 



a h--- a iJ°) ( i 

l<ii<-<i n <£ 

We define the conjugate vectors explicitly by the following: 

-i -l 



il+i2-\ \-in-nl+n(n— 1)/2 



(C.4) 



{£,n\ 



n 



q 



n{i-n) 



+-+in-ni+n(n-l)/2 



(C.5) 



Kii<-<»„<^ 



It is easy to show the normalization conditions [TS]: (£,n\ \ \ \£,n) = 1. In the massive regime 
where q = expr] with real r], conjugate vectors (£,n\\ are Hermitian conjugate to vectors \\£, n). 
Through the recursive construction fT27i2]) of P (£) s, it is easy to show the following [15] : 

P}Sj\£,n) = \\£,n), 



,n\\P, 



12— t 



(£,n\\. 



(C.6) 



Thus, the fusion projector is consistent with the spin-£/2 representation of U q (sl2}- 

In order to define Hermitian elementary matrices, we now introduce another set of dual 
basis vectors. For a given nonzero integer £ we define (£,n\\ for n — 0, 1, . . . , n, by 



n 



i 
n 



l<H<-<i„<^ 



+-+in)+n^-n(n-l)/2 



(C.7) 



They are conjugate to \ \£,n): (£, m\\ \ \£,n) = 5 mtn . 

In the massless regime where |g| = 1, matrix \\£, n)(£, n\\ is Hermitian: (\\£, n)(£, n\\y = 
\\£,n)(£,n\\. However, in order to define projection operators P such that PP = P, we define 
another set of vectors \ \£,n) in section 2.4. They are conjugate to the dual vectors (£, n\\. 
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D The massless fusion i?-matrices of the spin-1 case 



Let us evaluate the matrix elements of the massless monodromy matrix Tq^[ 2+ \\ ; £i), i.e. the 
massless L operator of the spin-1 representation. 



(2)0(1,1+), 



.(2)' 



5(2) 



;d.i) 



Here we have set inhomogeneous parameters = £1 and vffl = £1 — rj. Let us recall 



Rq U = Rq2^oi- For instance, we have Af 2 = A\ A\ + B^C^. In terms of b j = b(X — wf^) 



and c j = c(A — wf^) for j = 1, 2, the (1,1) element of A y ^' is given by 

(2, 1||A( 2+ )(A )ET) = (601 + &02 + <T 2 coiC 02 )/2 . 
Thus, setting w = A — £1, all the non-zero matrix elements of T^' 2+ \X ) are given by 



(2+) 



;d. 2 ) 



(2,0\\A {2 ^ 


(Ao)| 


|2,0) 


= (2,2||D( 2+ )(A )^2) = 1, 


(2,l\\A^ 


(Ao)j 


IM) 


= (2, 1||£>( 2+ >(A )||2, 1) = sinh(w + 77)/ sinh(w + 2?y) , 


(2,2||A (2V 


(Ao)j 




= (2, 0||£> (2+) (A )Eo) = sinhw/ sinh(w + 2rj) , 


(2,i||s( 2+ : 


(Ao)i 




= e~ u sinh 77/ sinh(w + 27?) , 


(2,2||B (2V 


(Ao)| 


|2,1> 


= [2] 9 g^ 1 e~ u sinh7//sinh(M + 2rj) , 


(2,0||C {2+ * 


(Ao)| 


|2,1> 


= [2] ? e" sinh 77/ sinh(-u + 2r\) , 


(2,1||C^ 


(Ao)| 


l?2) 


= qe u sinh 7// sinh(-u + 277) . 



(D.3) 

We should remark that the massive monodromy matrix T^[ 2+ \X ; £1) has the same matrix 
elements as the massless monodromy matrix T^[ 2+ \x ; £1). For instance, we calculate the (1,1) 
element of operator Af + ^ as follows. 



(2, lp (2+) (A )| |2, 1) = b m q~ 2 + 6 02 + <T 2 CoiC 2 = sinh(w + 77)/ sinh(w + 2rj) . 



;d.4) 



Let us define the matrix elements of the massless fusion R matrix of type (2, 2) as follows. 

^^(Ao-att = i(2,6i|U(2,6 ||<' 1 2+) (A -6)lKya|K)i. (D.5) 



Here ||2,c ) a and ||2,ci)i denote vectors in the auxiliary space and the quantum space 
V^ 2 \ respectively. Making use of matrix elements of the monodromy matrix of type (1, 2) we 
derive the fusion R matrix of type (2,2). For an illustration, let us calculate R^[ 2+ \u)li. 

a (2, 1||< 2+) (A - = (4 2+) (A)Cf +) (A - 77) + q- l c[ 2+ \X)Af + \x - 77)) q/[2] q . 

(D.6) 

(2+) (2+) (2) 

Evaluating operators A\ and C\ in the quantum space , we have 



(22+), uo _ [2] 9 e M sinh 7/ 
^o,i Woi- sinh( ^ + 2?7) - 



(D.7) 



The fusion i?-matrix becomes permutation n 0i i at u — 0. In fact, we have -Ro^^Woi = 1- 
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E Spin-s elementary matrices in global operators 



For integers ik and jk with 1 < i\ < ■ ■ ■ < i m < £ and 1 < ji < ■ ■ ■ < j n < £, we have 

-l 



\\£,m)(£,n 



II = 




' £ ' 




' £ ' 






m 


q 


n 



-(h+-+i m )+(ji+-+j n )+m(m+l)/2~n(n+l)/2 



?a(n^' n 4° II 

\k=l p=l;p^ix,... ,i m ,h>— Jn 9=1 



(E.l) 



Applying the spin-1/2 formulas of QISP [4J to (lE.li . we can express any given spin-s local 
operator in terms of the spin-1/2 global operators. It is parallel to the massive case P2]. Let 
us set %\ — 1, %i — 2, . . . , i m — m and j\ — l,j 2 — 2, . . . , j n = n in (lE.lj) . For m > n we have 



7? 



pg L n^ (i+) + D(i+) )K)ii D(i+) K 



i-l)^+fcj 



Q = l 



fc=l 



X 



] [ 5 (1+) K-i) 2s+fc ) [I I I (A^+D^)(w a ) P?} L . (E.2) 



fc=n+l 

For m < n we have 



fc=m+l 



•m, n (£ +) 



Hi 



n -1 



n 



(<-!)< 



a=l 



k=l 



x n c^ i+ \ w ^ l)2s+k ) n n (a^+d^k)^. (e.3) 

fc=m+l fc=m+l a=i£+l 

F Derivation of the density of string centers 

In terms of shifted rapidities A^ with A = 2s(a — 1) + a for a = 1, 2, . . . , N s /2 and a = 
1,2, ... ,2s, the Bethe ansatz equations for the homogeneous chain are given by 



M 



n ^a-xb+v) ^ forA=lj2j ... jM . 



sinh(AA + sr/) 

Smh(A\4 - ST)) J B =l;B^A sin H^A - A B - 77) 



(F.l) 



Putting A^ = fi a — (2s + 1 — 2a) and taking the product over a for a = 1, 2, ... , 2s, for the 
left-hand side of ( IF.lj) and for the right-hand side of (1F.1I) . we have 

Uivs inh ((fc- 1/2)77+^ 

/ yr I sinh(2sr/ - (/j a - // ft )) 2 -pr / sinhQr/ - (/i a - ^ b )) \ 2 . 

1 J 11 \ sinh(2 S r7 + (/i a -/i 6 )) 11 ^sinh^ + (// a - /i 6 )) ] f ~ 
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Taking the logarithm of (1F.2|) and making use of the following relation 

/ \ \ d 1 , /sinhf/cr? — A) \ .._ , 

we have the integral equation for the density of string centers, p(A), as follows. 

I oo / 28-1 \ 

p(A) = - J [ - ft) + J] 2K 2fc (A - /i 6 ) J P (A)rfA . (F.4) 



3=1 \ fc=i 

For 0<C<7r/mwe have the following Fourier transform: 



, smh(7ru;/2) 
Taking the Fourier transform of (IF .4ft we have the Fourier transform p(u;) of p(A) as follows. 

2s \ / 2s-l 



p(u) = J2K 2k ^(u) / 1 + ir 4s (^) + 2 #2*M 



1 



2cosh(Cw/2) 



v fe=i / \ fe=i 

Taking the inverse Fourier transform we obtain p(A) = l/2£ cosh(7rA/C). 

G Some formulas of the algebraic Bethe ansatz 

Applying the commutation relations between C and D operators we have 

M m M+l M+2 M+m 

(0\l[C(\ a )l[D(\ M+j ) = Yl E E G av .. am (X ir -- ,X M+m ) 

a=l j=l ai=l a,2=l;ai^ai a m =l;a\^a\,...,a m 

where 

n (\ \ \ T\(m\ s i n ngS^o^ sinh(A a . - A fe + V) \ 

G ai ---a m {Xl, - ■ ■ ,hl+ m ) = [[[diXa^iWjlL) — 2 — — — . (G.l) 

3=1 \ Llb=l;b^ ai ,..., aj smll lA aj . - X b ) J 

Let {\k}M be a set of Bethe roots. We have [H [5] 

M m M n M 

(01 n n ^k-) n s m i°>/<°i n n i°> 

k=l;k^ai,...,a m j=l 7=1 fe=l 7=1 

(-1 \ m 

(V({AJ)) *'({A a }\{A ai ,...,A am }U K,...,^}) J x J] 

x tt sinh(A a -Wj + 7]) -pj -pr sinh(A a - X a .) sinh(A aj . - A a J ^ 



46 



References 



[1] V.E. Korepin, N.M. Bogoliubov and A.G. Izergin, Quantum Inverse Scattering Method 
and Correlation Functions (Cambridge University Press, Cambridge, 1993) 

[2] N.A. Slavnov, Calculation of scalar products of wave functions and form factors in the 
framework of the algebraic Bethe ansatz, Theor. Math. Phys. 79 (1989) 502-508. 

[3] J.M. Maillet and J. Sanchez de Santos, Drinfel'd twists and algebraic Bethe ansatz, ed. M. 
Semenov-Tian-Shansky, Amer. Math. Soc. Transl. 201 Ser. 2, (Providence, R.I.: Ameri. 
Math. Soc, 2000) pp. 137-178. 

[4] N. Kitanine, J.M. Maillet and V. Terras, Form factors of the XXZ Heisenberg spin- 1/2 
finite chain, Nucl. Phys. B 554 [FS] (1999) 647-678 

[5] N. Kitanine, J.M. Maillet and V. Terras, Correlation functions of the XXZ Heisenberg 
spin-1/2 chain in a magnetic field, Nucl. Phys. B 567 [FS] (2000) 554-582. 

[6] N. Kitanine, J. M. Maillet, N. A. Slavnov and V. Terras, On the algebraic Bethe ansatz ap- 
proach to the correlation functions of the XXZ spin-1/2 Heisenberg chain, |hep-th/ 0505006, 

[7] M. Jimbo, K. Miki, T. Miwa and A. Nakayashiki, Correlation functions of the XXZ model 
for A < -1, Phys. Lett. A 168 (1992) 256-263. 

[8] M. Jimbo and T. Miwa, Algebraic Analysis of Solvable Lattice Models (AMS, Providence, 
RI, 1995). 

[9] M. Jimbo and T. Miwa, Quantum KZ equation with \q\ — 1 and correlation functions of 
the XXZ model in the gapless regime, J. Phys. A: Math. Gen. 29 (1996) 2923-2958. 

[10] T. Miwa and Y. Takeyama, Determinant Formula for the Solutions of the Quantum 
Knizhnik-Zamolodchikov Equation with |g| = 1, Contemporary Mathematics 248 (1999) 
377-393. 

[11] F. Gohmann, A. Kliimper and A. Seel, Integral representations for correlation functions 
of the XXZ chain at finite temperature, J. Phys. A: Math. Gen. 37 (2004) 7625-7651. 

[12] N. Kitanine, K.K. Kozlowski, J.M. Maillet, N.A. Slavnov, V. Terras, Algebraic Bethe 
ansatz approach to the asymptotic behavior of correlation functions, arXiv:0808.0227 
[math-ph] . 

[13] N. Kitanine, Correlation functions of the higher spin XXX chains, J. Phys. A: Math. Gen. 
34(2001) 8151-8169. 

[14] O.A. Castro-Alvaredo and J.M. Maillet, Form factors of integrable Heisenberg (higher) 
spin chains, J. Phys. A: Math. Theor. 40 (2007) 7451-7471. 



47 



[15] T. Deguchi and C. Matsui, Form factors of integrable higher-spin XXZ chains and the 
affine quantum-group symmetry, Nucl. Phys. B. 814 [FS] (2009) 405-438. 

[16] P.P. Kulish, N. Yu. Reshetikhin and E.K. Sklyanin, Yang-Baxter equation and represen- 
tation theory: I, Lett. Math. Phys. 5 (1981) 393-403. 

[17] H.M. Babujan, Exact solution of the one-dimensional isotropic Heisenberg chain with 
arbitrary spins s, Phys. Lett. A 90 (1982) 479-482. 

[18] H.M. Babujan, Exact solution of the isotropic Heisenberg chain with arbitrary spins: ther- 
modynamics of the model, Nucl. Phys. B 215 [FS7] (1983) 317-336. 

[19] A.B. Zamolodchikov and V.A. Fateev, A model factorized S'-matrix and an integrable 
spin-1 Heisenberg chain, Sov. J. Nucl. Phys. 32 (1980) 298-303. 

[20] K. Sogo, Y. Akutsu and T. Abe, New Factorized S-Matrix and Its Application to Exactly 
Solvable g-State Model. I; II, Prog. Theor. Phys. 70 (1983) 730-738; 739-746. 

[21] H.M. Babujan and A.M. Tsvelick, Heisenberg magnet with an arbitrary spin and 
anisotropic chiral field, Nucl. Phys. B 265 [FS15] (1986) 24-44. 

[22] A. N. Kirillov and N. Yu. Reshetikhin, Exact solution of the integrable XXZ Heisenberg 
model with arbitrary spin. I. the ground state and the excitation spectrum, J. Phys. A: 
Math. Gen. 20 (1987) 1565-1585. 

[23] T. Deguchi, M. Wadati and Y. Akutsu, Exactly Solvable Models and New Link Polyno- 
mials. V. Yang-Baxter Operator and Braid-Monoid Algebra, J. Phys. Soc. Jpn. 57 (1988) 
1905-1923. 

[24] M. Idzumi, Calculation of Correlation Functions of the Spin-1 XXZ Model by Vertex 
Operators, Thesis, University of Tokyo, Feb. 1993. 



[25] M. Idzumi, Correlation functions of the spin 1 analog of the XXZ model, hep-th/9307129 



[26] M. Idzumi, Level two irreducible representations of U q (sl2), vertex operators, and their 
correlations, Int. J. Mod. Phys. A 9 (1994) 4449-4484. 

[27] A.H. Bougourzi and R.A. Weston, Appoint correlation functions of the spin 1 XXZ model, 
Nucl. Phys. B 417 (1994) 439-462. 

[28] H. Konno, Free-field representation of the quantum affine algebra [/^(sLj) and form factors 
in the higher-spin XXZ model, Nucl. Phys. B 432 [FS] (1994) 457-486. 

[29] T. Kojima, H. Konno and R. Weston, The vertex-face correspondence and correlation 
functions of the eight-vertex model I: The general formalism, Nucl. Phys. B 720 [FS] 
(2005) 348-398. 



48 



[30] A. Caradoc, O. Foda and N. Kitanine, Higher spin models with domain wall boundary 
conditions, J. Stat. Mech. (2006) P03012. 

[31] L.A. Takhtajan, The picture of low-lying excitations in the isotropic Heisenberg chain of 
arbitrary spins, Phys. Lett. A 87 (1982) 479-482. 

[32] K. Sogo, Ground state and low-lying excitations in the Heisenberg XXZ chain of arbitrary 
spin S, Phys. Lett. A 104 (1984) 51-54. 

[33] F.C. Alcaraz and M.J. Martins, Conformal invariance and critical exponents of the 
Takhtajan-Babujian models, J. Phys. A: Math. Gen. 21 (1988) 4397-4413. 

[34] I. Affleck, D. Gepner, H.J. Schultz and T. Ziman, Critical behavior of spin-s Heisenberg 
antiferromagnetic chains: analytic and numerical results, J. Phys. A: Math. Gen. 22 (1989) 
511-529. 

[35] B.-D. Dorfel, Finite-size corrections for spin-S 1 Heisenberg chains and conformal properties, 
J. Phys. A: Math. Gen. 22 (1989) L657-L662. 

[36] L.V. Avdeev, The lowest excitations in the spin-s XXX magnet and conformal invariance, 
J. Phys. A: Math. Gen. 23 (1990) L485-L492. 

[37] F.C. Alcaraz and M.J. Martins, Conformal invariance and the operator content of the XXZ 
model with arbitrary spin, J. Phys. A: Math. Gen. 22 (1989) 1829-1858. 

[38] H. Frahm, N.-C. Yu and M. Fowler, The integrable XXZ Heisenberg model with arbi- 
trary spin: construction of the Hamiltonian, the ground-state configuration and conformal 
properties, Nucl. Phys. B 336 (1990) 396-434. 

[39] H. Frahm and N.-C. Yu, Finite-size effects in the XXZ Heisenberg model with arbitrary 
spin, J. Phys. A: Math. Gen. 23 (1990) 2115-2132. 

[40] H.J. de Vega and F. Woynarovich, Solution of the Bethe ansatz equations with complex 
roots for finite size: the spin S > 1 isotropic and anisotropic chains, J. Phys. A: Math. 
Gen. 23 (1990) 1613-1626. 

[41] A. Kliimper and M.T. Batchelor, An analytic treatment of finite-size corrections in the 
spin-1 antiferromagnetic XXZ chain, J. Phys. A: Math. Gen. 23 (1990) L189-L195. 

[42] A. Kliimper, M.T. Batchelor and P. A. Pearce, Central charge of the 6- and 19-vertex 
models with twisted boundary conditions, J. Phys. A: Math. Gen. 24 (1991) 3111-3133. 

[43] J. Suzuki, Spinons in magnetic chains of arbitrary spins at finite temperatures, J. Phys. 
A: Math. Gen. 32 (1999) 2341-2359. 



49 



[44] V.G. Drinfel'd, Quantum groups, Proc. ICM Berkeley 1986, pp. 798-820. 

[45] M. Jimbo, A g-Difference Analogue of U(g) and the Yang-Baxter Equation, Lett. Math. 
Phys. 10 (1985) 63-69. 

[46] M. Jimbo, A g-analogue of U(gl(N + 1)), Hecke algebra and the Yang-Baxter equation, 
Lett. Math. Phys. 11 (1986) 247-252. 

[47] M. Jimbo, Topics from representations of U q (g)- An Introductory Guide to Physicists, in 
Nankai Lectures on Mathematical Physics (World Scientific, Singapore, 1992) pp. 1-61. 

[48] J.M. Maillet and V. Terras, On the quantum inverse scattering problem, Nucl. Phys. B 
575 [FS] (2000) 627-644. 

[49] A. Kuniba, T. Nakanishi and J. Suzuki, Functional Relations in Solvable Lattice Models: I. 
Functional Relations and Representation Theory, Int. J. Mod. Phys. A 9 (1994) 5215-5266; 
II. Applications, Int. J. Mod. Phys. A 9 (1994) 5267-5312. 

[50] C.S. Melo and M.J. Martins, Algebraic Bethe Ansatz for £7(1) Invariant Integrable Models: 
The Method and General Results, Nucl. Phys. B 806 (2009) 567. 

[51] M.J. Martins and C.S. Melo, Algebraic Bethe Ansatz for £7(1) Invariant Integrable Models: 
Compact and non-Compact Applications. larXiv: 0902 .3476 . to appear in Nucl. Phys. B [FS] 
(2009) . 

[52] A. Nishino and T. Deguchi, The L(s^) symmetry of the Bazhanov-Stroganov model asso- 
ciated with the superintegrable chiral Potts model, Phys. Lett. A 356 (2006) 366-370. 

[53] M. Takahashi, Thermodynamics of One-Dimensional Solvable Models, (Cambridge Uni- 
versity Press, Cambridge, 1999). 

[54] F. Gohmann, N. P. Hasenclever and A. Seel, The finite temperature density matrix and 
two-point correlations in the antiferromagnetic XXZ chain, J. Stat. Mech. (2005) P10015. 



50 



